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ABSTRACT 

Chromium  doped  ^kermanite  crystals  were  investigated  by  optical 
spectroscopy  and  electron  paramagnetic  resonance  spectroscopy. 
The  results  were  compared  to  chromium  doped  forsterite  to  study 
the  roles  of  tetrahedrally  coordinated  chromium  ions  in  the  near- 
infrared  laser  transition. 

Akermanite  crystals  were  successfully  grown  by  Czochralski 
methods.  The  growth  processes  and  characteristics  of  various 
furnaces  were  explored  to  yield  good  optical  quality  crystals. 

The  energy  levels  of  chromium  ions  in  tetrahedral  sites  can  be 
explained  by  crystal  field  theory.  The  ground  state  energy  level 
splitting  under  an  external  magnetic  field  is  well  described  by 
spin-hamiltonian  matrices. 

From  electron  paramagnetic  resonance  measurements  of  the 
chromium  doped  ikermanite  and  forsterite  crystals  and  the 
similarity  of  their  optical  absorption  spectra,  tetravalent 
chromium  ions  were  believed  to  occupy  tetrahedral  sites  in  all 
the  samples.  However,  there  exist  considerable  differences  in 
zero-phonon  line  transitions  and  emission  spectra.  Therefore,  it 
has  not  been  completely  proved  that  the  mechanism  responsible 
for  near-infrared  laser  transition  in  chromium  doped  forsterite  is 
solely  due  to  chromium  ions  in  tetrahedral  sites. 
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Chapter  1 
Introduction 


1.1  Motivation 

Trivalent  chromium  ions  are  the  most  widely  used  dopant 
for  tunable  laser  crystals  such  as  the  chromium(Cr)  doped 
alexandrite  laser.  ^  A  chromium  activated  crystal,  forsterite 
(Cr3+;Mg2Si04),  was  recently  reported  to  lase  at  1250nm.2  Rager 
and  Weiser^  assigned  the  main  feature  of  the  absorption  spectra 
in  forsterite  to  Cr3+  ions  in  octahedral  sites.  Cr3+  ions  have  been 
found  by  electron  paramagnetic  resonance  methods  to  substitute 
for  the  two  nonequivalent  magnesium  octahedral  sites-the 
inversion  site  and  the  mirror  site.^ 

After  energy  level  calculations  were  performed,  it  was 
found  that  chromium  ions  in  octahedral  sites  do  not  have  the 
near-infrared  absorption  1 -1.2pm  which  was  observed  in 
Cr:Mg2Si04  .  Additionally,  the  large  absorption  bands  observed  in 
the  visible  frequency  regions  do  not  resemble  the  known 
absorption  bands  of  Cr^*^  in  other  octahedral  environments. 
Therefore,  it  is  unlikely  that  Cr3+  ions  in  octahedral  sites  are  the 
lasing  centers  in  Criforsterite.  Although  other  impurity  centers 
such  as  Fe3  +  ,  Fe^-*-  and  Cr2+  or  lattice  defects  such  as  Mg^  + 
vacancies  have  been  considered  to  be  the  lasing  centers,  they 
have  been  ruled  out  by  Verdun.5 


It  was  noted  that  the  near-infrared  absorption  strength 
scales  with  the  chromium  concentration  and  growth  atmospheres. 
This  scaling  indicates  that  chromium  is  responsible  for  the  laser 
transition;  however,  it  is  most  likely  some  valence  states  other 
than  Cr3  +  .  According  to  studies  by  Petricevic^  and  Verdun”^, 
crystals  grown  under  oxygen  atmospheres  have  absorption 
characteristics  which  can  be  explained  by  Cr^"*’  substituting  for 
silicon's  tetrahedral  sites  in  forsterite.  Time-resolved 
spectroscopy  indicated  the  existence  of  three  different  kinds  of 
emitting  centers;  which  can  be  attributed  to  Cr3+  and  Cr^'*’.®  All 
experimental  evidence  points  to  the  conclusion  that  ion  in 

tetrahedral  site  is  the  lasing  center. 

A  similar  situation  can  be  found  for  chromium  doped  garnet 
crystals  with  the  general  formula  A3B2C3O12  which  has  several 
substitution  sites  for  chromium.  The  trivalent  ion  is  in  the 

dodecahedral  site  with  coordination  eight,  while  is  in  the 

octahedral  site  with  coordination  six,  and  is  in  the  tetrahedral 

site  with  coordination  four.  Garnet  crystals  activated  with  rare 
earth  ions  such  as  Nd3+  doped  YsAlsOnfYAG)  are  important  solid 
state  laser  materials.  The  absorption  spectra  of  rare  earth  ions 
resembles  atomic  spectra,  with  many  narrow  lines  which  absorb 
little  flash-lamp  pumping  energy.  In  order  to  increase  the 
efficiency  of  the  Nd3+  laser  pumping,  chromium  ions,  with  their 
broad  absorption  bands,  are  co-doped  into  garnets.  Chromium 
normally  enters  the  garnets  as  Cr3+  in  the  octahedral  B  site  with 
an  absorption  bands  from  600  to  800nm.  But  it  was  noted  that 
there  is  also  an  additional  absorption  band  around  I -1.5pm  which 
could  not  be  attributed  to  Cr^'*'.  This  absorption,  which  overlaps 
the  lasing  range  of  Nd3+.  reduces  the  Nd^-*-  laser's  performance. 
This  absorption  was  initially  attributed  to  other  impurity  centers, 
but  was  finally  attributed  to  Cr^"*"  in  the  C  site  by  several 
authors. 9*^0  As  a  matter  of  fact,  this  impurity  center  has  already 
emerged  as  a  lasing  center  of  its  own.^  ^ 
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Due  to  the  presence  of  Cr  in  several  valence  states  as  well  as 
other  possible  types  of  charge  compensation  centers  in  the 
crystals,  there  will  be  overlapping  absorption  bands  from 
different  sites.  This  coupled  with  the  complex  structures  of  garnet 
crystals  and  forsterite,  makes  it  very  difficult  to  directly  identify 
the  impurity  center  responsible  for  the  near  infrared  absorption 
and  emission.  Although  the  center  has  been  assumed  to  be  Cr^'*’  in 
tetrahedral  sites  in  both  garnets  and  forsterite,  there  has  been  no 
clear  and  direct  evidence  of  Cr^'*'  ions  in  tetrahedral  sites. 

The  purpose  of  this  study  was  to  establish  the  existence  of 
chromium  ions  in  tetrahedral  sites  and  to  characterize  these 
centers. 


1.2  Identification  of  the  Chromium  Centers 

As  stated  in  the  last  section,  the  presence  of  octahedral  sites 
facilitates  the  existence  of  Cr3+  ions  in  forsterite.  This  will 
complicate  the  identification  processes.  In  order  to  isolate  the 
effects  solely  due  to  Cr*+  in  tetrahedral  sites,  another  kind  of 
crystal,  4kermanite  (Ca2MgSi207),  was  used.  Akermanite  is  a 
naturally  formed  mineral,  as  is  forsterite.  However,  it  is  different 
from  forsterite  in  that  there  are  no  octahedral  sites  in 
Akermanite.  The  chromium  dopant  can  only  enter  into  tetrahedral 
sites.  This  characteristic  makes  Akermanite  an  excellent  choice  to 
study  the  properties  of  chromium  ions  in  tetrahedral  sites. 

Even  with  Akermanite,  optical  methods  are  not  sufficient  to 
identify  the  impurity  center  due  to  the  broad  absorption  bands  of 
the  chromium  ions.  Therefore  electron  paramagnetic  resonance 
(also  called  electron  spin  resonance)  measurements  are  used  to 
assist  identification  of  the  center.  Electron  paramagnetic 
resonance(EPR)  can  directly  measure  the  ground  state  splitting  of 
the  dopant  ions,  which  can  yield  information  about  the  site 
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symmetry  and  the  valence  state  of  the  ion.  Based  on  the  data 
from  ikermanite,  chromium  doped  forsterite  can  be  analyzed 
more  easily  and  the  centers  can  be  readily  identified  by  the 
same  methods. 


1.3  Crystal  Growth  Facility 

In  order  to  grow  Ca2MgSi207,  a  radio  frequency(rO  induction 
heated  furnace  was  constructed.  The  alternating  current(AC) 
supplied  by  the  power  supply  produces  a  time-varying  magnetic 
field  which  in  turn  induces  an  AC  current  in  the  platinum  crucible 
wall.  According  to  Ohm’s  law,  the  crucible  will  be  heated  up. 

Owing  to  the  direct  coupling  between  the  crucible  and  the  rf 
coil,  the  furnace  has  a  large  temperature  gradient  and  is  suitable 
for  the  Czochralski  growth  method.  Before  growing  Ca2MgSi207. 
its  isomorphous  germanate,  Ba2MgGe207,  was  grown  due  to  its 
lower  growth  temperature  and  melt  viscosity.  Ba2MgGe207  was 
first  grown  in  our  lab  back  in  the  1960’s,  but  at  that  time  it  was 
used  as  a  laser  host  for  rare  earth  ions.  *  2 

After  Ba2MgGe207  was  successfully  grown,  Ca2MgSi207  was 
also  grown.  Due  to  some  difficulties  in  the  growth  of  Ca2MgSi207  in 
the  rf  induction  heated  furnace,  both  a  single  and  a  two-zone 
resistance  heated  furnaces  were  also  used  to  grow  Ca2MgSi207. 
These  furnaces  are  different  from  the  rf  furnace  in  that  they  are 
heated  up  by  thermal  convective  flow.  These  furnaces  have  a 
more  uniform  temperature  distribution  which  is  not  suitable  for 
Czochralski  growth;  however,  temperature  gradients  can  be 
obtained  by  adjusting  the  temperature  of  the  upper/lower  zones 
in  two-zone  furnace  or  by  using  an  air-cooled  seed  rod  in  the 
single-zone  furnace. 
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1.4  Contents  and  Organization  of  the  Thesis 

Chapter  two  will  describe  the  crystal  structure  of  both 
forsterite  and  4kermanite.  Their  structural  similarities  and 
differences  will  be  clarified.  Crystal  growth  processes  will  also  be 
explained  and  the  results  will  be  analyzed.  In  this  study,  three 
kinds  of  furnaces-  the  radio  frequency  induction  heated,  the 
single-zone  resistance  heated,  and  the  two-zone  resistance  heated 
furnaces  were  used.  Their  characteristics  are  explored  and 
compared. 

Chapter  three  establishes  the  theory  of  energy  levels  for 
chromium  ions  in  tetrahedral  sites.  It  is  mainly  based  on  crystal 
field  theory  which  treats  lattices  as  a  classical  distribution  of 
point  charges.  The  spin-hamiltonian  matrices  are  used  to  explain 
the  energy  level  splitting  of  the  ground  state,  which  is  essential  to 
interpret  the  electron  paramagnetic  resonance  spectra.  The  spin- 
orbit  interaction,  which  is  less  important  for  transition  metal  ions, 
is  also  generally  described. 

Chapter  four  lists  the  results  of  the  optical  spectroscopy  and 
electron  paramagnetic  resonance  spectroscopy.  Absorption  and 
emission  spectra  at  both  room  temperature  and  liquid  helium 
temperature  are  presented.  The  spectra  can  be  explained  by 
Tanabe-Sugano  diagram.  Also  the  EPR  field  for  resonance  spectra 
are  shown  and  can  be  explained  by  spin-hamiltonian  matrices  and 
the  space  group  symmetry  of  crystals. 

Chapter  five  summarizes  the  works  in  this  study  and  the 
comparisons  of  different  d^  configurations  in  different  host 
materials.  Finally,  possible  future  study  is  suggested. 
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Chapter  2 

Crystal  Preparation 


2.1  Crystal  Structure 

2.1.1  Forsterite  (Mg2Si04) 

Forsterite  belongs  to  the  olivine  group  with  orthorhombic 
symmetry  Pbnm.  The  unit  cell  contains  four  formula  units.  The 
unit  cell's  dimensions  are  a=4.7534A*  bsl0.1902A  and 

c=5.9783A.I  The  structure  contains  individual  [Si04]’^  tetrahedra 
linked  by  magnesium(Mg)  atoms  with  coordination  six.  The  same 
structure  with  different  divalent  atoms  linking  the  individual 
tetrahedra  constitute  a  series  of  compounds  (e.g.  fayalitefFeiSiOa), 
tephroite(Mn2Si04)  )  which  are  important  rock  forming  minerals 
of  the  earth. 

The  divalent  Mg  atoms  occupy  two  nonequivalent 
octahedral  sites  in  forsterite.  Half  of  the  Mg  atoms  are  located  at 
the  inversion  octahedral  site  with  symmetry  Ci,  the  other  half  at 
the  mirror  octahedral  sites  with  symmetry  Its  structure 
projected  on  (100)  plane  is  shown  in  Fig.  2.1.  The  crystallographic 
data  are  listed  in  Table  2.1. 

As  it  can  be  seen  from  Table  2.1,  the  silicon  tetrahedron  is 
distorted  such  that  it  has  only  mirror  symmetry. 
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Fig.  2.1  Unit  ceii  of  forsterite»  MgsSiOi.  Small  open  and 
solid  circles  are  Mg  atoms,  big  circles  are  O  atoms  and 
silicon  atoms  (not  shown)  are  at  the  centers  of  the 
tetrahedra.  (After  Ref.  2) 


Table  2.1  Crystallographic  data  of  forsterite.  (After  Ref.l) 


Z=4,  Symmetry  Pbnm(#62) 


Ion 

Site 

EHStM 

E1S9 

X 

y 

z 

KRB 

4a 

Q 

+2 

0 

0 

0 

ns 

4c 

Q 

+2 

0.99169 

0.27739 

0.25 

Si 

4c 

Q 

+4 

0.42645 

0.09403 

0.25 

Oi 

4c 

Q 

-2 

0.76549 

0.09156 

0.25 

02 

4c 

Q 

-2 

0.22164 

0.44705 

O} 

8d 

Cl 

-2 

0.27751 

0.16310 

0.03304 

The  Crrforsterite  crystals  used  in  this  study  were  provided 
by  Dr.  R.  Morris  of  Allied  Signal  Corporation.  These  crystals  were 
doped  with  a  sufficient  amount  of  chromium  to  distinctly  color  the 
crystals  olive  green.  When  looking  through  the  crystal’s  principal 
axes,  it  shows  greenish,  bluish  and  reddish  color  indicating  strong 
dichroic  effect. 


2.1.2  Akermanite  (CaiMgSiiO?  and  BaiMgCeiO?)^’** 

Akermanite  has  the  general  formula  A2BC2O7  and  belongs  to 
the  melilite  group  with  tetragonal  symmetry  P4  2im.  The  a2+  ions 
have  eight  fold  coordination,  forming  a  distorted  polyhedra  with 
two  different  A-0  distance.  The  B2+  ions  occupy  two  tetrahedral 
sites  (Tl)  located  at  the  corners.  ions  occupy  two  adjacent 

equivalent  tetrahedral  sites  (T2  and  T3)  to  form  [C207)6'  double 
pyramid  group.  So  the  crystal  structure  consists  of  [BOal^*  and 
[C20^]^•  parallel  to  the  [001]  direction  linked  by  A-0  bonds  as 
shown  in  Fig.  2.2.  The  structure  along  the  c  axis  is  shown  in  Fig. 
2.3. 


All  the  tetrahedral  sites  are  somewhat  distorted,  however. 
The  individual  symmetries  and  positions  for  Akermanite 
(Ca2MgSi207)  are  listed  in  Table  2.2.'^  All  the  atoms  are  arranged 
in  the  unit  cell  according  to  the  crystal  symmetry  P4  2im. 
Transition  metal  ions  are  suitable  to  substitute  for  the  smaller  B^*^ 
or  ions.  The  C^'*’  sites  which  we  are  interested  in  occupy  the 
4e  sites  in  the  notation  of  International  Tables  for  X-ray 
Crystallography.  Although  the  [COa]*’^  tetrahedra  are  distorted 
with  symmetry  Cs,  the  arrangements  of  the  4e  sites  make  them 
optically  equivalent  in  the  a-a  plane.  But  we  can  distinguish 
different  contributions  of  these  four  tetrahedra  in  EPR 
spectroscopy  by  applying  an  external  magnetic  field. 
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Table  2.2  Crystallographic  data  of  Ca2MgSi207.  (After  Ref. 
7) 


Z=2,  Symmetry  P4  2im(#113) 


Ion 

Site 

X 

y 

z 

Ca 

4e 

Q 

+2 

0.3355 

0.5-x 

0.5073 

Mr 

2a 

S4 

+2 

0 

0 

0 

Si 

4e 

Cs 

+4 

0.1396 

0.5-x 

0.9412 

Oi 

8f 

Q 

-2 

0.5 

0.5-x 

0.1804 

02 

4e 

-Q 

-2 

0.145 

0.5-x 

0.2583 

o, 

2c 

-2 

0.082 

0.182 

0.7909 

Akermanite-Ca2MgSi207  and  its  isomorphous  germanate 
Ba2MgSi207  were  grown  in  this  study.  The  tetrahedral 
coordinations  of  magnesium  are  unusual,  since  magnesium  in 
most  cases  is  found  in  coordination  six  such  as  in  the  case  of 
forsterite. 

Some  physical  properties  of  Ca2MgSi207  and  Ba2MgGe207  along 
with  forsterite  are  listed  in  Table  2.3. 


Table  2.3  Some  properties  of  &kermanite(Ca2MgSi207  and 

Ba2MgGe207)  and  forsterite.i*2,4,8,9 


composition 

unit  ceil 

refractive  index 

m.p. 

c/a/b(A) 

n(c)/n(a)/n(b) 

5.978/4.753/10.19 

1.670/1.635/1.651 

1890OC 

5.010/7.843 

1.640/1.632 

14540C 

5.56/8.36 

1.834/1.843 

~14030C 

The  local  symmetry  and  size  of  the  silicon  sites  of  forsterite 
and  Ca2MgSi207  are  similar.  The  germanium  site  in  Ba2MgGe207  is 
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a  little  larger.  The  interatomic  distances  and  site  symmetries  of 
these  two  compounds  are  shown  in  Fig.  2.4  and  listed  in  Table  2.4. 


Forsterite 

Mg2Si04 


a 


Average  Si-0 1.6346  A 


Akermanite 

Cd2MgSi20^ 

c 


Average  Si-0 1.633A 


Fig.  2.4  Interatomic  distances  and  site  symmetry  of  the 
silicon  site  in  Mg2Si04  and  Ca2MgSi207. 

Table  2.4  Interatomic  distances  of  silicon  tetrahedra.^’^^ 


Mg2Si04 

Ca2MgSi207 

Si-Oi 

1.6149A 

Si-Oi 

1.654A 

Si-02 

1.6530A 

Si-02 

I.592A 

Si-03 

1.6353A(twice) 

Si-03 

1 .643A(twice) 

O1-O2 

2.7434A 

O1-O2 

2.720A 

O1-O3 

2.7571A 

O1-O3 

2.562A 

O2-O3 

2.5527A 

O2-O3 

2.749A 

03-0’3 

2.5941A 

03-0'3 

2.6OOA 

Despite  its  similarity  with  forsterite,  the  distinctive  feature 
of  ikermanite  is  its  lack  of  octahedral  sites.  This  makes  it  an 
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excellent  material  to  study  the  effect  of  tetrahedral  coordination 
on  chromium  ions.  Additionally,  it  is  apparent  that  a  chromium 
ion  located  in  the  nearly  equivalent  average  Si-0  interatomic 
separations  (1.6346A  of  forsterite  vs  1.633A  of  Ca2MgSi20  7) 
should  lead  to  nearly  equivalent  crystal-field  strengths  and  thus 
transition  energies. 


2.2  Crystal  Growth 
2.2.1  Introduction 

The  phase  relations  of  the  Ca0-Mg0-Si02  ternary  system 
have  been  intensively  studied  because  numerous  compounds  are 
derived  from  this  system  (e.g.  forsterite(Mg2 S  i  O  4) . 
diopside(CaMgSi206),  merwinte(Ca3MgSi208)  and  etc.).  According 
to  the  phase  diagram  shown  in  Fig.  2.5,  Ca2MgSi207  is  congruently 
melting,  so  that  it  can  be  grown  from  a  charge  of  stoichiometric 
composition.  This  makes  it  a  good  candidate  for  Czochralski 
growth.  Its  corresponding  germanate  (Ba2MgGe207)  was  first 
grown  in  this  laboratory  in  the  I960*s  by  Czochralski  growth 
method. 8  This  crystal  also  congruently  melts  at  around  1403°C. 

Compared  to  Ba2MgGe207,  Ca2MgSi207  is  more  difficult  to 
grow  due  to  its  higher  melting  point;  in  addition,  molten  silicates 
have  high  viscosity  and  tend  to  solidify  as  glass.  In  considering 
these  reasons,  I  first  tried  to  grow  Ba2MgGe207,  then  Ca2MgSi207. 

300  gram  charges  of  powdered  BaC03  (Baker  Analyzed 
Reagent  grade),  MgO(Johnson  Matthey  Chemical)  and  Ge02(Eagle- 
Picher  Industrial)  were  mixed  together  according  to  their  molar 
prop  Mon  2:1:2.  4%  extra  germanium  oxide  was  added  in  order  to 
coir  ite  for  its  evaporation  during  the  growth  process. 
Chre  doping  was  provided  by  adding  Cr203  to  the  charge. 
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Fig.  2.5  Phase  diagram  of  system  CaO'SiOz* 
2Ca0  Mg0-2Si02.  (After  Ref.  11) 


Feed  materials  were  pre-heated  at  200°C  for  24  hours  to 
eliminate  the  water  adsorbed  by  these  materials.  Materials  were 
loaded  and  melted  into  a  250ml  platinum  crucible  repeatedly 
until  the  whole  300  gram  charge  could  fit  in.  The  pre-melting 
sequence  was 


200c 

-> 

3650C 

3 

hrs 

3650C 

-> 

9650C 

12 

hrs 

9650C 

-> 

lOOOoC 

1 

hrs 

lOOOoC 

-> 

12000c 

20 

hrs 

12000c 

-> 

1420OC 

11 

hrs 

14200  C 

-> 

200c 

10 

hrs 

Since  BaCOs  reacts  violently  to  release  CO2  in  the 
temperature  range  of  lOOO^C  to  1200oC,  the  step  must  be  very 
slow  in  this  range. 
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Both  a  induction  furnace  and  a  resistance  furnace  were  used 
to  grow  Ba2MgGe207.  The  characteristics  of  different  furnaces  will 
be  discussed  in  the  following  sections. 

The  basic  routine  for  a  Czochralski  growth  run  is^^ 

(1)  raise  the  temperature  a  few  degrees  higher  than  the  melting 
temperature. 

(2)  rotate  the  seed  rod  with  attached  seed  at  a  speed  of  15-60 
rpm  and  bring  slowly  into  contact  with  melt.  Then  adjust  the 
temperature  such  that  neither  growth  or  melting  of  the  seed 
occur. 

(3)  after  3-5  mins,  start  the  pulling  process  at  a  typical  rate  1-3 
mm/hrs. 

(4)  slowly  lower  temperature  to  achieve  a  desired  diameter. 

(5)  keep  adjusting  the  power  or  temperature  to  maintain  a 
constant  diameter. 

Using  the  above  processes,  undoped  and  0.01,  0.1  and  0.5 
mole  percent  chromium  (with  respect  to  Ge)  doped  Ba2MgGe207 
crystal  were  successfully  grown  with  diameter  1.5cm  and  length 
2-3  cm.  The  samples  were  confirmed  to  be  single  crystals  by 
observing  the  birefringence  of  the  crystals  under  polarizers.  The 
crystal  show  an  aqua-blue  color  and  have  a  weak  dichroic  effect. 
The  growth  direction  was  [001).  The  crystal  grown  by  using  the 
resistance  heated  furnace  shows  (100)  and  (110)  facet  planes.  On 
the  other  hand,  crystal  grown  by  the  rf  induction  heated  furnace 
has  no  facets,  which  is  probably  due  to  the  larger  radial  gradient 
in  the  rf  furnace.  After  X-ray  orientation,  the  boules  were  cut 
along  the  c  and  a  axes  to  obtain  samples  for  optical  measurements. 
To  facilitate  crystal  rotation  for  EPR  experiments,  crystals  were 
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cut  into  right  triangular  prisms  along  the  two  a-axes  and  the  c 
crystal  axis. 

After  growth,  the  crucible  can  be  cleaned  by  boiling  in  HCl. 

The  feed  materials  for  Ca2MgSi207  are  CaCOs,  MgO  and  Si02. 
They  are  mixed  also  according  to  the  molar  proportion  2:1:2.  The 
steps  of  preparation  of  Ca2MgSi207  samples  are  very  similar  to  the 
steps  of  Ba2MgGe207  except  minor  modifications  below: 

(1)  No  violent  chemical  action  when  CaCOs  releases  CO2,  so  the 
IOO-I2OO0C  heating  step  can  be  quicker. 

(2)  No  seed  was  available  at  the  beginning  of  growth  for 
Ca2MgSi207.  Seeds  were  obtained  from  material  nucleated  on  a 
platinum  wire.  This  process  can  be  done  by  quickly  immersing  the 
wire  in  the  melt;  otherwise,  instead  of  nucleation,  the  melt  will 
supercool  even  down  to  1 390^0 . 

(3)  Since  silicates  do  not  react  with  HCl,  the  best  way  to  clean  the 
crucible  is  to  dump  material  out  of  the  crucible  by  heating  the 
upside  down  crucible  to  the  melting  point  of  the  residue. 

By  using  this  standard  procedure,  undoped  and  doped  0.1 
and  0.25  mole  percent  chromium  (respective  to  silicon) 
Ca2MgSi207  were  successfully  grown  by  either  a  two-zone  or  a 
single  zone  resistance  heated  furnaces.  Cr:Ca2MgSi207  crystal 
shows  deep  blue  color  and  has  strong  dichroic  effect.  There  is  a 
strong  (001)  cleavage  plane  which  makes  seed  preparation  along 
[001]  direction  difficult.  So  most  of  growth  runs  were  performed 
along  either  [100]  or  [110]  direction.  The  resulting  crystals  shows 
strong  (100),  (110)  and  (001)  facet  planes  and  are  very 
unsymmetrical. 
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All  the  crystals  were  grown  in  ambient  atmosphere  in  order 
to  oxidize  Cr^-*-  to  Cr^*.  Typical  pulling  rates  range  from  0.5  to 
2mm/hrs  and  the  rotation  rate  between  15-60rpm. 

The  resulting  crystals  used  in  this  study  are  shown  in  Fig.  2.6. 


Fig.  2.6  Photo  of  Czchralski-grown  CrzCazM gS iiO 7, 
CriBazMgGeiO?  and  Cr:Mg2Si04  used  in  this  study. 

2.2.2  Crystal  Growth  by  a  RF  Induction  Heated  Furnace 

The  entire  system  includes  a  Pillar  Mark-6  lOKHz  30KW  rf 
power  supply,  a  Neslab  System  II  heat  exchanger  for  cooling,  a 
PC- AT  compatible  computer  with  appropriate  peripherals  for 
control  of  power/temperature,  pulling  and  rotation  speed.  The 
furnace  is  shown  in  Fig.  2.7. 


A.  Melt 

B.  Seed  and  Crystal 

C.  Platinum  Crucible 
0.  Alumina  Crucible 

E.  Alumina  Muffle 

F.  Alumina  Cap 


G.  RFColi 

H.  Thermocouple 

I.  Viewing  Window 

J.  ZIrcar  Insulation 

K.  Alumina  Cup 


Fig.  2.7  Schematic  of  the  rf  induction  heated  Czochralski 
furnace. 
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The  typical  axial  thermal  profile  of  this  rf  furnace  is  shown 
in  Fig.  2.8.  There  existed  a  large  thermal  gradient  of  about 
25oC/cm  above  the  melt  surface.  The  radial  gradient  was  also 
large.  Initially,  feedback  control  was  performed  through  readings 
of  a  thermocouple  near  the  melt  surface;  however  large 
temperature  variations  of  4-5oC  during  the  crystal  growth  run 
were  observed.  This  was  due  to  thermally  driven  convective  flow 
inside  the  furnace,  which  was  readily  seen  from  the  viewing 
window.  These  fluctuations  were  fed  back  to  the  controller  which 
in  turn  tried  to  correct  temperature  (by  correcting  the  power); 
which  caused  more  oscillations  in  the  temperature.  Instead  of 
temperature  control,  a  power  control  loop  was  used  while  an 
optical  pyrometer  was  used  to  measure  the  melt  temperature. 

The  power  control  loop  is  shown  in  Fig.  2.9.13  Power  can  be 
set  by  the  power  count(ranging  from  0  to  4095).  Via  a 
digital/analog  converter  and  an  optical  isolator,  the  power  count 
(PC)  is  converted  to  current  to  drive  the  power  supply.  Output 
voltage  of  a  power  measuring  bulb  (which  is  directly  proportional 
to  the  output  power)  is  read  and  digitized  for  the  computer.  This 
measurement  is  used  to  correct  the  power  count  to  get  the  a 
constant  power  output. 

The  algorithm  for  correcting  the  power  count  is  shown  as 
followed. 


(PC)next=(PC)iast*(  1+P8*(  l'P®'*'®*’mcasurcd/powerscipoint))  (2.1) 

The  parameter  Pg  (power  gain)  was  chosen  with  respect  to 
the  measurement  response  time  and  the  given  sampling  interval. 
The  sampling  interval  is  1  sec  in  this  system.  Although  this  is  not 
a  linear  control  system,  it  works  like  a  PI  controller  with  steady 
state  error  equal  to  zero.  After  several  tries,  Pg=0.2-0.3  gave 
satisfactory  results  with  power  regulation  to  within  0.1%.  The 
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temperature  (°C) 


temperature  reading  on  top  of  the  melt  continued  to  fluctuate 
about  2-3  degrees.  The  fluctuation  is  believed  to  result  from  the 
thermally  driven  air  flo*^.  Fig,  2.10  and  Fig.  2.11  show  the 
corresponding  power  and  temperature  fluctuations  for  the 
temperature  control  and  the  power  control  mode.  The  advantage 
of  using  power  control  rather  than  temperature  control  can  easily 
be  seen. 

By  using  the  control  methods  described  above,  several 
0.01%Cr:Ba2MgGe207  crystal  boules  were  pulled.  During  growth,  it 
was  observed  that  it  was  very  difficult  to  manually  control  the 
diameter  of  the  crystal.  When  expansion  of  the  crystal  diameter 
was  visually  noticed  during  the  temperature  cooling  and  pulling 
was  started,  the  crystal  emerging  from  the  melt  already  had  a 
large  diameter.  It  appears  that  the  crystal  had  started  growing 
inside  the  melt.  Additionally,  most  crystals  grown  showed  ring 
shaped  bubble  inclusions  along  the  boules  in  the  regions  between 
the  center  and  the  edge.  This  behavior  can  be  explained  by  a  non¬ 
flat  growth  interface.  It  is  well  known  that  a  proper  balance  of 
forced  convection  (due  to  seed  rod  rotation)  and  natural 
convection  (due  to  temperature  gradient  of  the  melt)  are  crucial  to 
produce  the  desired  growth  interface,  and  good  optical  quality 
crystals.  Further  work  on  the  flow  pattern  of  the  melt  needs  to 
be  done  to  achieve  good  crystal  quality. 


2.2.3  Crystal  Growth  by  a  Two-Zone  Resistance  Heated 
Furnace 

Due  to  the  difficulty  of  controlling  a  rf  furnace  with  the 
previous  setup,  a  resistance  heated  furnace  seemed  to  be  an 
attractive  alternative.  Thus,  a  two-zone  furnace  shown  in  Fig.  2.12 
was  used  to  grow  undoped  and  chromium  doped  Ca2MgSi207.  The 
control  system  is  very  similar  to  that  of  the  rf  furnace;  except  that 
the  control  thermocouple  is  located  at  the  center  of  the  lower 
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Fig.  2.10  Measured  power  and  temperature  fluctuation 
the  temperature  controi  mode. 
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Fig.  2.12  Schematic  of  the  two-zone  resistance  hi 
furnace. 
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heating  element.  The  thermal  gradient  can  be  set  by  adjusting  the 
relative  power  input  to  the  upper  and  lower  heating  elements. 

By  controlling  the  temperature  of  the  lower  heater 
thermocouple(LHT),  a  temperature  variation  of  ±0.0100  can  be 
obtained.  The  same  O.OloC  variation  was  observed  on  both  the 
upper  heater  thermocouple  and  a  thermocouple  located  above  the 
melt;  which  illustrates  the  more  steady  thermal  environment  in 
this  kind  of  furnace. 

By  using  the  technique  described  in  section  2.2.1,  undoped 
seed  material  (shown  in  Fig.  2.13)  was  nucleated  on  a  platinum 
wire.  Several  similar  foot-shaped  crystals  were  pulled  later.  But 
generally  speaking,  the  growth  was  not  very  successful.  These 
crystals  either  grew  rampantly  or  melted  all  the  way,  showing 
that  the  temperature  near  the  growth  interface  was  either  too 
high  or  too  low.  This  unusual  behavior  puzzled  me  for  quite  a 
while.  Originally,  I  thought  that  seed  was  reacting  with  the 
aiumina(Al20 3)  seed  rod.  But  when  I  put  a  platinum  button 
between  the  seed  and  the  seed  rod,  the  same  phenomenon 
occurred.  Therefore,  this  is  not  the  case.  A  detailed  axial 
temperature  measurement  was  performed  by  setting  LHT=1596oC 
both  with  and  without  a  crucible  to  determine  possible  causes. 
The  results  are  shown  in  Fig.  2.14  and  Fig.  2.15.  Although  the 
temperature  measurement  contained  a  lot  of  noise,  it  can  give  us 
a  general  trend  of  the  temperature  variations  along  the  axis. 

When  the  crucible  was  not  in  the  furnace,  the  temperature 
decreased  monotonously  from  the  bottom  to  the  top.  However, 
when  the  crucible  was  put  in,  there  existed  an  inverted 
temperature  gradient  in  the  first  1cm  above  the  melt  surface.  Also 
in  the  latter  case,  the  furnace  needed  less  power  input  to  achieve 
the  same  LHT  temperature.  The  facts  stated  above  show  that  with 
the  crucible  in  place,  the  outside  wall  of  the  crucible  can  reflect 
heat  back  to  the  heating  elements  and  somehow  heat 
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Fig.  2.13  Photo  of  undoped  CaiMgSiiO?  grown  from  a 
platinum  wire. 


distance  from  melt  surface  (mm) 

Fig.  2.14  Temperature  profile  along  furnace  axis  without 
a  crucible. 
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Fig.  2.15  Temperature  profile  along  furnace  axis  with  a 
crucible. 


can  transfer  from  the  seed  rod  to  the  seed.  Under  these  conditions, 
a  good  seeding  temperature  was  obviously  difficult  to  achieve. 

2.2.4  Crystal  Growth  by  a  Single*Zone  Resistance  Heated 

Furnace 

A  single  zone  resistance  heated  furnace  (shown  in  Fig.  2.16) 
contrary  to  that  of  a  rf  induction  heated  furnace  has  a  rather 
uniform  temperature  distribution  inside  the  furnace.  In  fact,  due 
to  the  heat  conduction  path  through  the  bottom,  temperature 
decreased  from  the  melt  surface  to  the  bottom  of  the  crucible.  A 
small  inverted  temperature  gradient  also  existed  above  the  melt 
surface,  most  likely  for  the  same  reason  as  in  the  two-zone 
furnace. 
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.16  Schematic  of  the  singie*zone  r< 
furnace. 
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This  isothermal  condition  is  not  suitable  for  Czochralski  growth.  In 
order  to  produce  a  temperature  gradient  above  the  melt  surface, 
an  air-cooled  seed  rod  was  used  to  force  a  desired  gradient.  After 
air  cooling  was  provided,  controlled  growth  was  achieved. 

From  one  of  the  samples  grown,  a  possible  vapor  etching  of 
the  seed  was  observed.  In  order  to  solve  this  problem,  two  tricks 
were  used: 

(1)  before  growing,  use  the  cooling  rod  to  blow  air  above  melt 
surface  for  1-3  mins. 

(2)  prior  to  pulling,  dip  the  seed  into  the  melt  to  fix  the  broken 
regions  inside  the  grooves  which  are  used  to  attach  the  seed  to 
the  platinum  seed  holder. 


CajMgSijO^:  Cr 

4-:0-9l 
10.6SU 

l_J _ I _ I _ I _ I _ ^ 

0  I  2  3  4  5  cm 

MIT 

Lab  for  Advanced  Solid  State  Laser  Mati 
_ _ _ 


Fig.  2.17  Photo  of  the  crystal  grown  by  using  air  blowing 
and  crystal  dipping  techniques. 
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After  these  two  tricks  were  performed,  crystal  growth  proceeded 
rather  smoothly  and  yielded  the  crystal  shown  in  Fig.  2.17.  This 
crystal  was  grown  rather  slowly  at  about  0.8mm/hrs  compared  to 
2mm/hrs  in  the  rf  furnace’s.  Some  defects  were  incorporated  in 
the  middle  region  and  unsymmetrical  growth  was  observed 
indicating  that  the  crucible  was  not  properly  centered.  Seed 
orientation  along  [100]  or  [110]  induced  faceted  growth. 

2.2.5  Conclusion 

(1)  The  characteristics  of  three  different  kinds  of  furnaces  were 
explored.  The  rf  induction  heated  furnace  has  a  greater  thermal 
gradient  and  is  best  controlled  by  power  input.  The  resistance 
heated  furnace  need  external  cooling  to  force  a  temperature 
gradient  near  the  growth  interface. 

(2)  All  the  crystals  were  grown  from  a  crucible  less  than  half-full 
of  melt.  This  low  height  of  melt,  which  inhibits  both  natural  and 
forced  convective  flows,  will  produce  an  undesired  growth 
interface  and  thus  poor  optical  quality  crystals. 

(3)  Growing  along  [001]  and  a  properly  centered  crucible  are 
needed  to  attain  symmetric  growth. 

(4)  Defect  inclusions  made  crystals  unacceptable.  Some  growth 
parameters  need  to  be  changed  such  as  rotation  rates, 
temperature  lowering  rates,  etc.. 


40 


References 


1  K.  Fujino,  S.  Sasaki,  Y.  Takeuchi  and  R.  Sadanaga,  "X-ray 
Determination  of  Electron  Distributions  in  Forsterite,  Fayalite 
and  Tephroite,"  Acta  Cryst.,  B37,  p513,  1981. 

2  W.  A.  Deer,  R.  A.  Howie,  J.  Zussman,  "Rock  Forming  Minerals," 
VI,  P2-33,  John  Wiley  and  Sons  Inc.,  New  York,  1962. 

3  J.  V.  Smith,  "Reexamination  of  the  Crystal  Structure  of 
Melilite,"  American  Mineralogist,  V39,  p643,  1953. 

4  W.  A.  Deer,  R.  A.  Howie,  J.  Zussman,  "Rock  Forming  Minerals," 
VI,  P236-255,  John  Wiley  and  Sons  Inc.,  New  York,  1962. 

5  S.  F.  Bartram," Crystal  Structure  of  Y2SiBe207,"  Acta  Cryst., 
B25,  p791,  1969. 

6  T,  H.  Allik,  M.  J.  Ferry,  R.  J.  Reeves,  R.  C.  Powell,  W.  W.  Hovis, 
D.  P.  Caffey,  R.  A.  Utano,  L.  Merkle  and  C.  F.  Campana, 
"Crystallography,  Spectroscopic  Analysis,  and  Lasing 
Properties  of  Nd3+:Ba2ZnGe207."  J.  Opt.  Soc.  Am.,  B7(7),  p 
1190,  1990. 

7  Structure  Reports,  V17,  p574,  1953. 

8  M.  Alam,  K.  H.  Gooen,  B.  D.  Barlolo,  A.  Linz,  E.  Sha^,  L. 
Gillespie  and  G.  Jannery,  "Optical  Spectra  and  Laser  Action  of 
Neodymium  in  a  Crystal  Ba2MgGe207",  J.  Appl.  Phys.,  V39, 
p4728,  1968. 

9  N.  A.  Sirazhiddinov,  N.  N.  Mirababaeva,  R.  G. 
Grebenshchikov,  and  E.  V.  Strganov,  "Isomorphism  of 
Barium  Germanosilicates  in  the  Ba2Si04,  Ba2Ge04llMSi03, 
MGeOa  Quaternacy  Reciprocal  System,"  Russian  Journal  of 
Inorganic  Chemistry,  VI 9(6),  p817,  1974. 

10  J.  R.  Smith,  R.  M.  Hazen,  "The  Crystal  Structures  of  Forsterite 
and  Hortonolite  at  Several  Temperatures  Up  to  900oC," 
Amer.  Mineral.,  V58,  p588,  1973. 


41 


11  J.  B.  Ferguson  and  H.  E.  Merwin,  "The  Ternary  System  CaO- 
Mg0-Si02,"  Amer.  Joum.  Sci.,  4th  ser.,  V48,  pl09,  1919. 

12  J.  C.  Brice,  "Crystal  Growth  Processes",  Chapter  4,  Blackie  & 
Son  Ltd.,  Bishopbriggs,  1986. 

1  3  Peter  H.  Kloumann,  "Operations  Manual  for  PC  Based  Crystal 
Growth  System,"  Crystal  Physics  Lab.,  1991. 


42 


Chapter  3 

Theory  of  Chromium  Ions  in 
Tetrahedral  Sites 

3.1  Environment  Consideration 

When  a  small  amount  of  chromium  metal  ions  are 
introduced  into  solids,  we  can  assume  that  the  chromium  ions  only 
interact  with  the  local  environment  (  in  this  case,  oxygen  ions). 
The  main  hamiltonian  for  the  impurity  ion  is 

"  1  ,  Ze*  "  e* 

«=  I{— Pk - +  )*i,  •  Sk  +  2:  — +  V(r  )}  (3.1) 

k=i  2m  k<x 

The  first  term  is  the  electron  kinetic  energy  and  the  second  term 
represents  the  Coulomb  interaction  due  to  electrons  and  the 
central  nucleus;  the  third,  forth  and  fifth  terms  respectively  are 
the  spin-orbit  interaction,  electron-electron  interaction  and 
interaction  due  to  electrons  and  local  environments.  The  first  four 
terms  have  the  same  form  as  for  the  free  ion,  but  the  last  term  is 
solely  due  to  the  interaction  of  the  ion  with  the  rest  of  solid.  In 
order  to  write  down  the  explicit  form  of  the  fifth  term,  we  must 
make  several  assumptions. 
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First,  we  can  model  the  environment  as  a  classical 
distribution  of  charge.  Second,  the  ionic  crystal  model  which 
assumes  that  the  electrons  are  solely  bounded  to  metal  ions  or  the 
cations  can  be  used.  Thus  we  can  write  down 


Vc,ysta|('‘K)  =  ?7® 

I  ll 


kl 


(3.2) 


which  we  call  the  crystal  field. 

The  relative  magnitudes  of  the  five  terms  in  Eqs  3.1  are 
important  for  calculating  the  energy  levels,  especially  the  spin- 
orbit  interaction  and  the  crystal  field  interactions.  According  to 
thek  relative  strengths,  they  can  be  classified  the  following  two 
ways. 

(1)  weak  field  scheme:  The  crystal  field  is  smaller  than  the  spin- 
orbit  interaction.  We  therefore  solve  for  the  free  ions'  energy 
levels  first  and  then  treat  the  crystal  field  as  a  first  order 
perturbation.  Rare  earth  or  actinide  ions  are  usually  in  this  class, 
because  4f  and  5f  electrons  are  well  shielded  from  the  crystal 
field  by  other  electrons. 

(2)  strong  field  scheme:  The  crystal  field  is  stronger  than  the  spin- 
orbit  interaction.  So  we  must  consider  the  crystal  field  first  and 
then  treat  spin-orbit  interaction  as  a  first  order  perturbation. 
Transition  metal  ions  in  oxide  crystals  are  usually  in  this  class. 
The  d  electron  orbitals  extend  to  the  exterior  of  the  atoms  and 
have  strong  interaction  with  the  crystal  field. 

We  would  expect  that  chromium  ions  in  oxides'  tetrahedral 
sites  should  belong  to  the  latter  case.  From  now  on  only  the 
strong  field  scheme  will  be  considered. 
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In  the  strong  field  scheme,  we  first  solve  the  crystal 
potential  component  and  then  the  Coulomb  repulsion  between  the 
electrons  (if  the  metal  ion  has  more  than  one  electron  which  is  the 
case  for  divalent,  trivalent  and  the  tetravalent  chromium).  Finally, 
spin-orbit  coupling  and  low  symmetry  components  of  the  crystal 
field  are  treated  as  perturbations. 

The  crystal  potential  possess  (and  also  the  electron -electron 
interaction  and  kinetic  energy)  all  the  symmetry  elements  of  the 
group  of  the  site.  For  a  tetrahedral  site,  the  appropriate  group  is 
Td.  Before  chromium  ions  are  put  into  the  solid,  the  free  ion  has 
full  rotational  symmetry.  The  characters  of  the  angular 
momentum  state  t  under  rotation  a  degrees  about  the  z  axis  and 
rotation  a  plus  inversion  are^ 


sin(/-Kl/2)a 

X  (a)  = - 

sln(1  /  2)a 

/  sin(£  + 1  /  2)a 

(la)  =  H)  - 

8in(1  /  2)a 


(3.3) 


After  the  ions  are  put  into  tetrahedral  sites  (shown  in  Fig.  3.1), 
the  full  rotation  symmetry  of  free  ions  is  reduced  to  the 
symmetry  of  Td,  whose  operations  consist  of: 

•  identity  E 

•  8  C3  about  4  body  diagonals  corresponding  to  rotation  ±27c/3 

•  3  C2  about  X,  y,  z  directions 

•  6  S4  about  X,  y,  z  corresponding  to  rotation  of  in/2 

•  6  Od  planes  that  are  diagonal  reflection  planes 


The  character  table  for  Td  symmetry  is  shown  in  Table  3.1. 


Fig.  3.1  Symmetry  operations  of  tetrahedral  sites. 

Table  3.1  Character  table  for  Td  symmetry  and  the 
corresponding  fs2  full  rotation  state. 
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Full  ROT 

for  /=2 
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We  can  also  use  Eqs.  3.3  to  calculate  the  corresponding 
characters  for  the  free  orhit  angular  momentum  state  t-1  (d  state) 
under  the  operations  of  id  symmetry.  The  results  are  shown  in 
Table  3.1. 

According  to  the  reduction  theorem  of  reducible 
representations,  the  characters  of  reducible  representation  x(Ck) 
(full  rotation  in  tetrahedral  sites)  is  a  linear  combination  of 
characters  for  the  irreducible  representation  of  group  Td. 

Comparing  the  characters  of  full  rotation  in  a  tetrahedral 
site  with  the  characters  of  Td  group,  we  get 

X(/=2)  =  E  +  T2 

In  Td  symmetry,  the  five  fold  degenerate  d  functions  will 
therefore  split  into  two  sets  with  irreducible  representation  E  and 
T2.  And  the  corresponding  eigenfunctions  for  the  H  crystal-field 
are 


E  dx2.y2=(  1  /^/2)(d2+d.2)=(^/3/2)(x2-y2) 

dz2=do=l/2(3z2-r2)=l/2(2z2-x2-y2) 

T2  dxy=(l/V2)(d2-d.2)=>/3xy 

(l/V2)(di-d.i)=%/3xz 

dyz=(l/V2)(di+d.i)=V3yz  (3.4) 


Although  group  theory  describes  the  degeneracy  and 
splitting,  it  does  not  tell  us  the  magnitude  and  the  orderings  of  the 
splitting.  These  quantities  can  be  calculated  from  the 
eigenfunctions  and  crystal  fields.  The  potential  of  identical 


nearest  neighbors  for  various  arrangements  of  configurations 
are^ 


V 


bcc 


=  -  4  —  Y,,  cJlZ.rs  +  Vsz^  ) 


40 


40 


44 


for  8  ligands  in  the  vertices  of  a  cube  and  the  ion  sitting 

in  the  body  center 

“'•I'oo  -2— 


tet 


40 


40 


44 


for  4  ligands  in  the  vertices  of  a  regular  tetrahedron 


VT 


V 


oct 


=  ♦m's: 


40 


40 


c 

"44' 


for  6  ligands  in  the  vertices  of  a  regular  octahedron 


(3.5) 


where 

4  =  1/2)  I 

_  4ne  - 

"  2k+lJ  rI'+I 

p  ( R)  is  the  spatial  charge  distribution 


After  calculating  the  hamiltonian  matrices  associated  with  the 
basis  functions  E  and  T2  using  the  tetrahedral  potential,  we  get 
the  following  results: 


<T2IV,e,IT2>=l/V45(4yQQ+(8/9)Y4  q?) 
<EIVtt,IE>=l/^/4i(4Ygo-(4/3)Y4 

<T2lV„tlE>=0  (3.6) 


48 


which  confirm  the  energy  splitting  results  derived  from  the  group 
theory. 

The  difference  of  energy  between  the  two  representations  is 

E(T2)-E(E)=20/(9•^/4S)Y4Q?=10Dq  (3.7) 

Dq  is  an  important  parameter  in  the  crystal  field  theory  which 
describe  the  magnitude  of  the  cubic  crystal  field.  And  from  Eqs. 
3.5,  it  is  very  easy  to  get  the  following  relations  between 
different  cubic  environments  assuming  that  the  charge 
distributions  around  the  ions  and  the  distances  between  the  ions 
are  the  same. 


(Dq)bcc=2(Dq)tct='2(4/9)(Dq)oct  (3.8) 


3.2  and  d^  in  Cubic  Environments 

If  there  exists  more  than  one  d  electron  in  the  atom,  we 
must  consider  the  Coulomb  interaction  between  d  electrons.  But 
first,  we  can  use  the  direct  product  of  groups  to  get  the  irreducible 
representations  of  many-clectron  wavefunctions.  For  example, 
Cr^+  have  two  d  electrons.  If  both  electrons  are  in  the  E  level. 
Coulomb  interaction  will  split  it  into  Ai,  A2  and  E  levels  by  the 
following  formulas. 


E  X  E  =  Aj+A2  +  E 


By  the  same  reason. 


For  spin  part. 


'^2  ^  ^2  ~  ^1  ^  ^1  ^2 
E  X  T2  =  Tj  + 
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(S=l/2)  X  (S=l/2)  =  S=0  or  I 


(3.9) 


then  we  can  apply  the  well-known  formula  to  construct  the 
wavefunctions: 

0(SM)-Znii,m2  0(simi)0(s2m2)<simiS2ni2lSM>  (3.10) 

<sim  iS2ni2lSM>  are  Wigner  coefficients  (or  Clebsch-Gordan 
coefficients)  which  are  tabulated  in  standard  group  theory  books. 
After  applying  Eqs.  3.10,  some  of  the  combinations  vanish  and 
only  the  following  terms  survive. 

1Ai(T22),  1Ai(E2),  1E(T22),  1E(E2).  1T2(T22),  1T2(T2E),  3Ti(T22), 
3Ti(T2E).  1Ti(T2E),  3T2(T2E)  and  3A2(E2) 

The  number  in  the  upper  left  hand  corner  represents  the  spin 
degeneracy  2S+1,  while  the  representation  inside  the  bracket  is 
the  original  electron  configurations  without  Coulomb  interactions. 

Since  the  Coulomb  interaction  shown  in  Eqs  3.1  also  has  the 
full  symmetry  of  the  crystal,  energy  levels  with  the  same 
symmetry  will  mix  together.  Of  ail  the  possible  cross  terms 
between  levels,  we  only  need  to  consider  terms  having  the  same 
symmetry.  The  appropriate  matrix  elements  have  been  compiled 
by  Tanabe  and  Sugano.3  These  matrix  elements  and  eigenenergies 
are  listed  in  the  Appendix.  According  to  the  eigenenergies 
calculated  from  the  formula  in  the  Appendix,  we  can  construct 
the  Tanabe-Sugano  (T-S)  diagrams  that  represent  the  energy  ratio 
E/B  as  a  function  of  Dq/B  (  the  definition  of  B  is  in  Appendix).  The 
Tanabe-Sugano  diagrams  which  have  achieved  great  success  in 
interpreting  the  spectra  of  transition-metal  ions  in  cubic 
environments  will  be  used  in  conjunction  with  the  absorption 
spectra  in  Chapter  3.  A  typical  energy  levels  ordering  for  Cr^"*" 
ions  in  tetrahedral  sites  are  shown  in  Fig.  3.2. 
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For  cl3  configuration,  we  can  follow  exactly  the  same 
procedure  by  using  the  direct  product  method.  The  resulting 
symmetry,  matrix  elements  and  eigenenergies  are  also  listed  in 
Appendix. 

It  is  interesting  compare  the  d**  with  dlO‘“  configuration.  A 
dlO-n  configuration  forms  a  full  shell  in  conjunction  with  a  d^^ 
configuration.  As  a  matter  of  fact,  the  states 
t26*n(S  iri)e^‘™(S2r2)Sr-Mr  are  complementary  states  to  the 
t2"(Siri)e®(S2r2)SrMr.  Under  the  approximation  that  all  E  and  T2 
functions  are  d-functions,  the  matrix  elements  in  N-6,7,8,9  can  be 
obtained  from  N=l,2,3,4  by  replacing  t2‘^(Siri)e“(S2r2)SrMr  with 
t2^‘“^S iri)e'^*"'(S2r2)SrMr.  So  in  the  d*  configuration,  we  can 
still  use  the  matrix  elements  and  eigenenergies  of  the  d^ 
configuration. 


Fig.  3.2  Energy  level  diagram  for  a  d^  ions  in  a 
tetrahedral  field. 
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3.3  Spin*Orbit  Interaction  and  Low  Symmetry  Fields 

The  energy  splittings  due  to  the  spin-orbit  interaction  and 
the  low  symmetry  fields  in  transition  metal  ions  are  usually  one 
order  of  magnitude  smaller  than  the  electronic  Coulomb 
interaction.  Therefore,  we  can  treat  the  spin-orbit  interaction  and 
the  low  symmetry  field  as  perturbations. 

The  hamiltonian  associated  with  the  spin-orbit  interaction 
can  be  written  as^ 


H-Si^(ri)li*si  (3.11) 

which  after  some  arithmetic  manipulations  can  be  rewritten  by 
using  tensor  operator  equivalence  as 

K,o-l/V2[-V+ia(lT|)+iV.,p(lTi)l+ 

l/V2t-V.a('Ti)+iV.ip(lTi)l+V0Y(>Ti)  (3.12) 

where  VtiaaZjSiillia,  Vtip*LiSjtllip.  Vjor*'2iSi0tiy 

Si±l=+1  /^^(Sixtisiy),  SiO=Siz 
tia=^(ri)Iix»  ti3=^(ri)Iiy,  tiY=^(ri)liz 

One  way  to  look  at  the  Eqs.  3.12  is  that  the  orbital  angular 
momentum  states,  lix.  liy,  liz.  transform  as  the  Ti  irreducible 
representation  in  cubic  symmetry.  Therefore,  we  can  use  the  basis 
functions  of  Ti  to  evaluate  the  hamiltonian. 

In  real  crystals,  there  is  always  some  distortion  to  perfect 
cubic  environments,  so  we  often  get  a  low  symmetry  component 
of  the  crystal  field.  For  a  low  symmetry  crystal  field,  we  can  also 
write  down  the  corresponding  tensor  operators  based  on  the  form 
of  potential.^  For  example,  in  Dan  symmetry,  the  potential  can  be 
expressed  as  Vu(Eg),  and  the  D3  d  symmetry  as 
{V^(T2g)+VT,(T2g)+V(;(T2g)} 
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Thus  we  can  calculate  the  matrix  elements  between 
different  terms  as  <asr  i  H $ o+H|ow-symmctry  I  a  '  s  '  r  '  >  .  Then, 
standard  perturbation  procedure  can  be  used  to  get  the  new 
energy  levels. 

Sometimes,  the  low  symmetry  field  will  have  larger  strength 
than  the  Coulomb  interaction;  in  these  cases  the  perturbation 
procedure  cannot  be  used.  Under  these  circumstances,  we  must 
first  calculate  the  energy  splitting  due  to  the  whole  crystal  field, 
then  calculate  the  Coulomb  energy,  and  lastly  treat  spin-orbit 
interactions  as  perturbations.  The  whole  procedure  is  very  similar 
to  the  procedure  stated  in  section  3.1,  but  the  lower  symmetry  is, 
the  more  complex  the  procedure  is.  Otsuka  has  worked  out  the 
case  of  tetragonal  fields. ^ 


3.3.1  Energy  Level  Splitting  of  Ground  State 

The  method  stated  in  the  last  section,  although 
straightforward,  is  still  very  complicated,  especially  when  the 
spin-orbit  interaction  and  low  symmetry  field  parameters  are 
included  in  the  calculations.  So  instead  of  solving  the  real 
hamiltonian,  we  can  approach  this  problem  phenomenologically 
by  constructing  an  effective  hamiltonian. 

In  this  section,  I  will  use  the  spin  hamiltonian  matrices  to 
account  for  the  effects  of  all  the  low  order  interactions,  which 
include  the  spin-orbit  interaction  and  spin-spin  interaction^  on 
the  d^  ground  state  energy  of  chromium  ions  in  tetrahedral  sites, 
(this  calculation  will  be  used  in  electron  spin  resonance 
experiments  in  Chapter  4).  For  Cr^'*'  ions  in  forsterite  and 
4kermanite,  the  ground  state  is  (see  the  Appendix). 

Due  to  the  presence  of  the  low  symmetry  field,  the  ground 
state  is  usually  an  orbital  singlet  with  representation  or 
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2s+1A2.  The  ground  state  wavefunctions  will  be  split  further  by 
the  spin-orbit  interaction,  spin-spin  interaction,  etc..  The  general 
form  of  the  ground  state  hamiitonian  for  s=l/2  is 


Kso+Iow  symm  = 


^  a 

^c+id 


c~id^ 


which  can  be  rewritten  by  using  the  well-known  spin  matrices 
0(l/2,±l/2)(notice, 

J^so+low  symm  =l/2(a+b)l+(a-b)Sz+2cSx+2dSy  (3.13) 


ri  0^ 

fo 

ro 

n 

where  1  = 

Sy  = 

oj 

Sz=  „ 

-V 

lo  V 

Oy 

J 

lo 

In  the  same  manner,  we  can  express  the  S=1  system  by  linear 
combinations  of  Syy  SySz^Sz^y?  ^z Sx^^x^z> 

SxSy+SySx>  Sz**  Sx^'Sy*.  Initially,  it  might  seem  that  we  do  not 
get  any  advantage  as  we  still  have  a  lot  of  linear  coefficients  to 
determine.  But  by  noticing  the  following  facts:  (1)  the  orbit 
singlet  ground  state,  'F  (^5  + ^  A  i  ,2M  e  i  ,2)  has  the  same 
transformation  properties  as  ©(SM)f(Ai  or  A2);  (2)  the  orbital 
part  transforms  as  the  symmetry  of  the  wavefunctions(Ai  or  A2); 
(3)  then  the  base  functions  of  the  spin  matrix  @(SM)  must 
transform  as  Ai  (the  identity  irreducible  representation),  in  order 
for  the  whole  wavefunctions  to  have  the  symmetry  Ai  or  A2;  (4) 
the  hamiitonian  has  time-reversal  symmetry^,  we  can  greatly 
simplify  the  effective  hamiitonian. 

For  example,  for  a  spin-triplet  state  in  D2  symmetry 
environment,  we  can  use  0(1M)  (1,  Sx>  Sy,  Sz,  SySz+SzSy, 
SzSx+SxSz»  SxSy+SySx»  Sz**  Sx*-Sy*)  as  the  bases.  The  time 
reversal  operator  is  defined  as  Ks0(l/2,m)=(l/2)-'n0(l/2,-m).  By 
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applying  the  time  reversal  symmetry  to  0(1M),  the  Sx,  Sy,  Sz 
linear  terms  which  change  sign  are  required  by  symmetry  to 
have  zero  coefficients.  Also  by  applying  the  2-fold  rotational 
symmetry  along  x,  y,  z  axes,  the  terms  SySz+SzSy,  SzSx+SxSz, 
SxSy+SySx  which  change  sign  should  also  have  zero  coefficients. 
Then  the  effective  hamiltonian  can  be  simplified  as 

H  =  D[Sz"-l/3S2]+E(Sx*-Sy2)  (3.14) 

where  D  and  E  are  the  fine  structure  parameters  which  need  to  be 
determined  experimentally.  The  constant  term  -1/3S2  is  included 
so  that  the  trace  of  the  matrix  is  zero. 

For  states  with  symmetry  T  other  than  Ai  and  A2,  we  must 
use  the  bases  (pCT'Y)  which  transform  like  the  direct  product  of 
r  X  r.  For  example,  the  hamiltonian  of  excited  state  3X2  can  b® 
written  as^ 


Jt=-A(Lz^*2/3)-r(Lx^*Ly3)+(X.xSxLx+^ySyLy+  XzSzLz)  (3.15) 

Let  us  go  back  to  Eqs.  3.14.  The  whole  hamiltonian  is  (where  I 
neglect  the  common  factor  8^/3) 

(DI3  0  E  ^ 

0  -20/ 3  0 

,  E  0  D/3, 

After  solving  the  above  matrix,  the  eigenenergies  are  found  to  be 

Ei=D/3+E,  E2=D/3-E,  E3=-2D/3  (3.16) 

As  can  be  seen,  the  ground  state  will  split  into  triplet  states  with 
the  energies  given  above. 
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Now,  consider  the  case  when  an  external  magnetic  field  is 
applied  to  the  crystal  (this  is  the  case  in  electron  paramagnetic 
resonance  experiments).  According  to  quantum  mechanics,  the 
magnetic  moments  produced  by  the  electron  spin  and  orbital 
motion  will  precess  around  the  magnetic  field  H  with  the  energy 

AE  =  -p-H  =  .pHcos(p,H)  (3.17) 

where  p  is  magnetic  moments  due  to  the  spin  and  orbital 
magnetic  moments 

li=-  (PB/ft)(2iLi+22iSi)=-  (pb/»)(L’+2S')  (3.18) 

where  pb  Is  Bohr  magnetron,  and  L’,  S',  the  total  orbit  and  spin 
angular  momentum. 

But  when  the  ground  state  is  an  orbital  singlet  state  due  the 
low  symmetry  field,  its  orbital  angular  momentum  is  zero.  The 
reason  is: 

The  wavefunction  of  nondegenerate  ground  state  must  have  the 
form  of  Eqs.  3.4  to  match  the  corresponding  irreducible 
representation,  (notice  the  wavefunctions  of  Eqs.  3.4  which 
transform  as  xy,  yz,  zx,  Sr^-x^-y^,  x2-y2  are  linear  combinations  of 
angular  momentum  eigenstates,  d.2,  d-i,  do,  di  and  d2).  So  the 
wavefunctions  have  real  forms.  But  the  operator  Lz  for  the  z 
component  of  the  orbital  angular  momentum  is  a  purely 
imaginary  operator. 


L2=:ift(x8/ dy-yd/ 3x)  (3.19) 

Lz  is  also  a  hermitian  operator  whose  eigenvalue  is  real.  So  the 
only  possible  solution  of  the  eigenequation  is  fz=0.  In 

this  case,  the  orbital  angular  momentum  is  quenched.  So  in  the 
ground  state,  we  can  ignore  the  spin-orbit  interaction  and  also  let 


L-0  in  Eqs.  3.17.  So  by  combining  Eqs.  3.17  and  3.14,  we  get  the 
spin-hamiltonian  for  a  ground  state  of  3a  j  or  3a2  in  an  external 
magnetic  field. 

K  =  pH.g.S+D[SzJ-l/3S2]+E(Sx*-Sy*)  (3.20) 

where  p=p.B/ft  and  the  g  tensor  is  used  instead  of  the  isotropic  g 
value  of  2. 

The  energies  of  the  three  states  as  a  function  of  the  external 
magnetic  field  are  plotted  in  Fig.  3.3.  From  the  Fig,  3.3,  we  can 
notice  that  the  magnitude  and  the  ordering  of  the  splitting  will  be 
different  depending  on  the  direction  of  the  magnetic  field.  This 
special  feature  is  extremely  useful  to  identify  the  site  symmetry 
of  the  dopant  ions  in  the  solid.  Detail  will  be  discussed  in  section 
4.3. 


Fig.  3.3  ^Ai/^A.2  level  splitting  as  a  function  of  H 
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Chapter  4 

Spectroscopy  and  Analysis 

4.1  Experiment  Procedure 

4.1.1  Optical  Spectroscopy 

Polarized  absorption  spectra  can  be  obtained  by  using  an 
available  Perkin-Elmer  Lambda  9  spectrophotometer  in 
conjunction  with  polarizing  prisms.  Emission  spectra  are  recorded 
by  using  a  McPherson  218  Scanning  Monochromator  with  786A 
Digital  Scanning  Control  and  an  EG&G  5301  Lock-In  Amplifier. 
Excitation  of  the  emission  is  by  one  of  the  following  sources:(l) 
Spectra  Physics  171  Argon  Laser;  (2)  argon  laser  pumped 
Ti^+rsapphire  laser;  (3)  argon  laser  pumped  dye  laser;  (4)  Spectra 
Diode  Lab  SDL-2422  GaAlAs  semiconductor  laser.  The  excitation 
source  is  chosen  depending  on  the  absorption  band  pumped.  The 
whole  set  up.  is  shown  in  Fig.  4.1. 

For  fluorescence  decay  or  lifetime  measurements,  (1)  a 
Holobeam  QG  500  pulsed  frequency  doubled  Nd3+:YAG  laser;  (2)  a 
pulsed  NdzYAG  pumped  Tirsapphire  laser;  (3)  a  pulsed  Molectron 
DL-II  series  nitrogen  laser  pumped  dye  laser  (also  depend  on  the 
absorption  band  pumped)  were  used  for  excitation.  The  decaying 
signal  is  recorded  by  an  Analogic  Data  6000  Universal  Waveform 
Analyzer  in  conjunction  with  a  Tektronix  7623  oscilloscope.  The 
whole  set  up  is  shown  in  Fig  4.2. 
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Fig.  4.1  Experimental  set-up  for  emission  measurements. 


Fig.  4.2  Experimental  set-up  for  life-time  measurements. 
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Detectors  ranging  from  IR/visible  photomultipliers,  Ge,  Si, 
and  InSb  photodetectors  are  used  depending  on  the  sensitivity 
required  and  the  spectral  range  of  interest. 

In  order  to  resolve  some  fine  structures,  low  temperature 
spectroscopy  is  needed.  An  Air  Products  and  a  Janis  Cyrogenic 
variable  temperature  cryostats  are  used  in  conjunction  with  the 
Perkin-Elmer  Lambda  9  spectrophotometer  and  emission  set-up 
to  obtain  liquid  nitrogen  or  helium  temperature  spectra. 


4.1.2  Electron  Paramagnetic  Resonance  Spectroscopy 

Due  to  the  broad  energy  bands  of  transition  metal  ions  in 
solids  (which  result  from  lattice  interactions),  it  is  difficult  to  use 
optical  spectroscopy  to  identify  the  impurity  center.  On  the  other 
hand,  the  electron  paramagnetic  resonance  (EPR)  spectroscopy  has 
the  following  advantages^; 

(1)  The  microwave  frequency  used  in  EPR  are  10^  to  10^^  Hz 
which  correspond  to  10*^  to  10*^  eV.  Such  a  small  energy  enables 
us  to  investigate  the  small  splitting  due  to  the  low  symmetry 
field  and  the  external  magnetic  field  which  is  not  accessible  by 
optical  methods. 

(2)  By  changing  the  direction  of  magnetic  field,  different  splittings 
due  to  the  local  symmetry  of  the  impurity  center  can  be  explored. 

(3)  By  changing  the  strength  of  magnetic  field,  different 
magnitude  of  energy  splitting  can  be  produced. 

(4)  Due  to  different  nuclear  magnetic  dipole  moments  of  different 
atoms,  sometimes  we  can  identify  the  impurity  by  observing  the 
characteristics  of  its  hyper  fine  structure. 
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(5)  By  observing  the  temperature  dependence  of  resonance  and 
using  the  pulsed  EPR  method,  spin-spin  and  spin-lattice 
relaxations  of  the  impurity  can  be  studied. 

Because  of  the  above  advantages  and  also  the  lack  of  an 
effective  way  to  identify  the  impurity  center  by  optical  methods, 
EPR  spectroscopy  of  Cr:forsterite  and  Cr:&kermanite  crystals  were 
performed  to  assist  the  identification,  (this  part  of  the 
experimental  work  was  undertaken  by  graduate  student  M.  H. 
Whitmore  and  Prof.  D.  J.  Singel,  Chemistry  Department,  Harvard 
University). 2 

The  procedure  of  EPR  experiments  can  be  found  in  any 
standard  magnetic  resonance  textbooks.^*'^ 


4.2  Results  and  Analysis  of  Optical  Spectroscopy 
4.2.1  Cr:Forsterite 

The  room  temperature  polarized  absorption  spectra  of 
Cr;Mg2Si04  (grown  under  an  oxidizing  atmosphere)  are  shown  in 
Fig.  4.3.  The  three  distinct  absorption  peaks  at  571,  662,  and 
742nm  are  found  and  also  weaker  absorption  bands  in  the  near- 
infrared  region  and  at  400nm.  Since  all  these  bands  have 
considerable  overlap,  the  interpretation  of  these  spectra  is 
difficult.  By  observing  the  absorption  strength  vs.  growth 
atmospheres,  Petricevic^  and  Verdun^,  Baryshevski7  et  al,  have 
assigned  the  various  bands  to  Cr3+  in  octahedral  sites  and  Cr^'*’  in 
tetrahedral  sites. 

In  my  sample  (grown  under  an  oxidizing  atmosphere),  the 
main  features  correspond  to  those  identified  for  Cr^‘*‘  ions  in  silicon 
sites.  Before  looking  into  the  spectra,  it  is  worthwhile  to 
understand  the  local  environment  of  the  Si  site  in  forsterite.  The 
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Fig.  43  Room  temperature  absorption  spectra  of  Q’:Mg2SiQt. 
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Fig.  4.4  Symmetry  descending  and  associated  selection  rules. 

Si  site  as  shown  in  Table  2.1  has  symmetry  Cs  lowered  from  the  Td 
symmetry  of  a  regular  tetrahedron.  The  distorted  tetrahedron 
was  shown  in  Fig.  2.4.  As  a  first  order  approximation,  we  can  treat 
the  tetrahedron  having  strain  along  a  axis  which  lower  the 
symmetry  to  Csv  The  progression  of  descending  symmetry(T(i- 
>C3v->Cs)  and  selection  rules  for  electric  dipole  transitions  can  be 
derived  from  group  theory^  and  are  shown  in  Fig.  4.4.  From  the 
figure,  we  notice  that  3A2'>^Ti  is  the  only  allowed  transition  in 
the  Td  symmetry.  Therefore,  we  can  expect  that  this  3a2->3Ti 
transition  will  be  the  strongest  for  Cr^'*'.  The  three  large  peaks  in 
the  visible  region  should  therefore  belong  to  the  3a2*>^Ti 
transition.  The  much  weaker  near>infrared  absorption  is  due  to 
3A2->^T2,  as  can  be  seen  from  the  energy  diagram.  By  using  the 
mean  value  of  the  three  peaks  and  the  cubic  field  approximation 
(see  the  Appendix),  the  crystal  field  parameters  are  determined 
to  be  Dq=-1010.8cm'i,  B=569.7cm-l.  The  parameter  C  is  set  to  be 
4.4B  as  in  the  free  ion  case  due  to  the  lack  of  reliable  transition 
assignments  to  the  E  states.  By  using  these  values  for  Dq,  B,  C,  and 
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also  the  eigenenergy  formula  in  the  Appendix,  the  energy  levels 
of  Cr4+  in  forsterite  were  calculated.  The  corresponding  energy 
diagram  and  Tanabe-Sugano  diagram  are  shown  in  Fig.  4.5. 

Let  us  compare  the  result  with  the  data  published  by 
Petricevic^,  Verdun^  and  Moncorge^.  The  Dq's  values  are  very 
similar.  However,  there  exist  considerable  differences  for  the  B 
and  C  values.  This  is  due  to  the  different  assignments  of  the 
centers  of  absorption  bands  to  their  associated  irreducible 
representations.  Even  at  low  temperature,  absorption  bands  still 
overlap  such  that  the  assignments  of  absorption  bands  are 
extremely  difficult.  But  by  noticing  that  the  electric-dipole 
selection  rule  for  S  is  AS=0,  we  expect  that  the  dominant 
transitions  should  be  3A2->3T2  and  ^Ti.  From  the  known  Ni2'‘’:MgF2 
spectra  (which  is  a  d^  configuration  in  octahedral  sites;  and  the 
energy  structure  of  d^  in  octahedral  sites  should  be  very  similar  to 
that  of  d2  in  tetrahedral  sites^),  the  absorption  around  350nm, 
374nm  and  444nm  should  result  from  transition  from  3A2  to  the 
second  ^Ti.  Even  so  the  more  precise  calculation  is  very  difficult 
for  the  following  reasons:  (1)  the  lack  of  reliable  data  from  3A2  to 
a  spin-singlet  state;  (2)  the  large  low-symmetry  field  (can  be  seen 
from  the  large  splitting  of  ^Ti  and  3T2  states);  (3)  absorption 
bands  overlapping  with  Cr3+  absorption  bands;  and  (4)  low  local 
symmetry  Cs  of  the  silicon  site.  Better  energy  level  calculations 
have  been  obtained  either  by  the  tetragonal  field 
approximation^.lO  or  by  treating  low  symmetry  fields  as 
perturbations.il  But  neither  method  has  considered  the  following 
internal  limitations  of  the  crystal  field  theory: 

(1)  in  a  non  inversion  symmetry  environment,  the  configuration 
mixing  with  p  states  and  even  f  states  may  be  large  enough  that 
we  can  not  ignore  it. 
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(2)  in  such  a  high  valence  state,  we  would  expect  that  some  part 
of  the  bonding  between  Cr-0  will  be  covalent.  This  contradicts  the 
basic  ionic  crystal  approximation  of  the  crystal  field  theory.  In 
fact,  silicon  atoms  in  forsterite  have  electronic  charges  +2.1 
determined  from  effective  distribution  radii  ^2  instead  of  the 
formal  charge  4+. 

For  the  illustrative  purpose,  the  cubic  field  approximation  is 
good  enough  to  explain  the  general  characteristics. 

For  E//a,  we  can  clearly  observe  the  zero-phonon  line  and 
the  corresponding  phonon  sidebands  at  17521cm*l(570nm), 
18244cm-H548nm)  and  18980cm*U527nm).  The  average  phonon 
energy  for  this  structure  is  730cm*  K  In  order  to  observe  the  fine 
structure  of  other  polarizations  or  frequency  regions,  liquid- 
helium  temperature  absorption  experiments  were  performed. 
The  lOK  absorption  spectrum  for  E//b  (the  same  axis  as  for  the 
dominant  emission)  is  shown  in  Fig.  4.6.  Again,  three  peaks  at 
13474cm*l(742nm),  I4217cm*l(704nm)  and  14935cm*H668nm) 

are  observed.  The  average  phonon  energy  is  730cm*  ^  and  is 
consistent  with  the  E//a  room  temperature  absorption  spectrum. 
The  zero  phonon  line  for  the  3a2->^T2  transition  is  also  easily 
identified  at  1090nm  with  very  clear  vibronic  structures.  This 
result  is  very  similar  to  Ni2  +  (d8)  in  MgF2.  Therefore,  the 
assumption  that  Cr^'*'  ions  are  in  silicon  sites  is  consistent. 

Forsterite  has  emission  in  the  near-infrared  region.  The 
room  temperature  polarized  emission  spectra  for  the  near- 
infrared  region  excited  by  0.3W  argon  laser  pumping  at  488nm 
with  k//c  are  shown  in  Fig.  4.7.  The  spectra  were  recorded  using  a 
McPherson  0.3  monochromator  with  a  1200  grooves/mm  grating 
and  a  Ge  detector.  For  E//a,  the  small  emission  signal  may  be  due 
to  the  slight  mis-alignment  of  the  crystal.  For  the  E//c  spectrum, 
we  notice  that  there  is  an  emission  region  extending  below 
lOOOnm  which  is  quite  different  from  the  E//b  spectrum.  For  the 


68 


OO 

r:’- A 

>-  II  = 
UHUI 


(  U13)  *udoo  uojidjosqv 


Fig.  4.6  lOK  absorption  spectrum  of  Cr:Mg2St04  for  E//b. 
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Fig.  4.7  Room  temperature  polarized  emission  spectra  of 
Cr:Mg2Si04. 
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emission  longer  than  llOOnm,  the  fluorescence  decay  is 
exponential  with  a  lifetime  2,7ps  at  room  temperature.  But  for  the 
emission  below  llOOnm,  a  double  exponential  decay  was 
observed  indicating  that  two  kinds  of  emission  centers  exist.  The 
additional  emission  center  is  attributed  to  Cr3+  in  Mg(Cs)  site.7  The 
decay  curve  at  1200nm  is  shown  in  Fig.  4.8. 

There  is  little  doubt  that  the  near-infrared  emission  with  a 
lifetime  2.7ps  and  near-infrared  absorption  belong  to  the  same 
electronic  transition,  especially  after  we  compare  the  low 
temperature  emission  and  absorption  spectra  in  the  NIR  region. 
They  show  an  almost  perfect  mirror  image  with  the  pure 
electronic  transition  located  at  I090nm  as  shown  in  Fig.  4.9.  The 
strong  phonon  assisted  transition  peaks  in  the  emission  spectra 
show  that  at  least  two  different  low  frequency  phonons  with 
energy  -150  cm*^  and  490cm*  1  are  involved,  although  the  exact 
assignment  is  not  very  clear.^*^^ 


time(^s) 

Fig.  4.8  Cr:Mg2Si04  room  temperature  fluorescence  decay 
at  1200nm. 
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Fig.  4.9  Near  infrared  E//b  absorption  and  emission  of 
Cr:forsterite  at  liquid  helium  temperature. 
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4.2.2  Cr;Ca2MgSi207 


The  room  temperature  polarized  absorption  of 
0.25%Cr:Ca2MgSi2O7  is  shown  in  Fig.  4.10.  Since  Ca2MgSi207  has 
tetragonal  symmetry  (compared  to  orthorhombic  symmetry  of 
forsterite),  only  the  E//c  and  E±c  spectra  are  distinct.  When 
compared  to  the  spectra  of  forsterite,  we  can  immediately  identify 
the  similar  absorption  peaks  at  579nm,  660  nm  and  767nm  and 
also  the  weaker  absorption  bands  in  the  NIR  region. 

But  the  absorption  bands  around  440  and  660nm  which  we 
attribute  to  Cr^-*-  in  octahedral  sites  are  absent  in  Cr:Ca2MgSi207 
crystal.  This  result  is  consistent  with  the  crystal  structure  of 
Ca2MgSi207,  because  as  stated  in  Section  2.1,  Ca2MgSi207  does  not 
have  any  octahedral  sites  and  also  Cr^'*’  ions  should  enter  the  Si 
sites  where  is  no  need  for  any  charge  compensation. 

The  lOK  absorption  spectra  are  shown  in  Fig.  4.11.  Based  on 
the  spectra  and  the  same  calculation  procedure  as  in  the  last 
section,  the  cubic  field  parameters  can  be  obtained  by  using  the 
formula  in  Appendix  to  be  Dq=-996cm*l,  B=546.5cm*^  and  C  was 
assigned  to  4.4B.  The  calculated  and  observed  spectra  are  shown 
in  Fig.  4.12.  Although  the  detailed  splitting  can  not  be  determined 
by  this  procedure,  the  energy  diagrams  reasonably  describe  the 
general  characteristics  of  Cr^'*'  in  the  Si  site  of  Ca2MgSi207. 

The  low  temperature  absorption  spectra  do  not  show  any 
structure  at  all  either  in  the  visible  or  NIR  frequency  region 
except  an  additional  peak  at  555nm(180S0cm'i).  This  should  be  a 
phonon  replica  for  the  peak  at  579(  1 7270cm*  0.  showing  the 
phonon  energy  780cm* ^ 

Two  other  samples  with  concentration  1  and  0.1  mole 
percent  with  respect  to  the  Si  were  also  examined.  All  samples 
have  the  same  spectral  characteristics. 
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Fig.  4.10  Room  temperature  polarized  absorption  spectra 
of  Cr:Ca2MgSi207. 
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Fig.  4.  11  lOK  polarized  absorption  spectra 
CriCazMgSiiO?. 
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Fig.  4.12  Calculated  and  observed  energy  levels  of  Cr^  * 
ions  in  Ca2MgSi207. 

The  effort  of  finding  emission  from  Cr:Ca2MgSi207  was 
unsuccessful.  The  crystal  was  pumped  along  the  a  axis  with 
sample  length  0.715cm  by  using  an  Ar+  laser  at  488nm.  Several 
other  pumping  sources  such  as  a  dye  laser  at  650nm  (with 
Rhodamine  6G  dye)  and  a  TiiSapphire  Laser  ranging  from  700- 
850nm  were  also  used.  Still  no  emission  was  detected.  The 
measurement  was  attempted  at  liquid  nitrogen/helium 
temperature  but  no  signal  was  detected.  Based  on  the  sensitivity 
of  the  spectrometer,  the  Cr:Ca2MgSi207  emission  intensity  must  be 
at  least  a  factor  of  100  less  than  that  of  Cr:forsterite. 

In  order  to  check  the  possible  influence  of  a  valency  change 
in  the  Ca2MgSi207  crystals  on  the  spectra,  I  attempted  to  reduce 
the  crystals  by  annealing  in  an  atmosphere  with  an  oxygen  partial 
pressure  of  10"  atm  at  750^0  for  36  Hours.  The  room 
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temperature  absorption  remained  the  same  and  emission  was  still 
absent  (indicating  low  oxygen  diffusion).  Lower  oxygen  partial 
pressures  were  obtained  with  forming  gas  (95%  N2  and  5%  H2)  and 
the  samples  were  heated  to  1000°  C  for  24  hours  and  then 
quenched.  After  reduction,  the  absorption  of  Cr:Ca2MgSi207  was 
smaller  by  -25%  throughout  the  spectra  indicating  that  some  Cr4  + 
may  have  reduced  to  Cr3+.  But  no  new  features  appeared,  which 
needs  to  be  investigated  further.  This  Cr^+  absorption  variations 
with  annealing  atmospheres  are  similar  to  Cr4+  in  forsterite.  After 
annealing,  Cr3+  ions  in  forsterite  may  further  reduced  to  Cr^'*’  to 
give  extra  absorption  bands  as  observed, 


4.2.3  CrtBaiMgGeiO? 

Fig.  4.13  shows  the  polarized  absorption  spectra  obtained  for 
0.5%Cr:Ba2MgGe2O7.  Their  spectral  shape  is  similar  to  that  of 
Cr:Ca2MgSi207  and  Cr:forsterite,  although  red-shifted  several  tens 
of  nm.  The  spectra  consist  of  rather  large  peaks  at  636,  728  and 
832nm  and  a  weaker  NIR  absorption. 

When  the  absorption  spectra  was  measured  at  liquid  helium 
temperature  as  shown  in  Fig.  4.14,  several  structures  showed  up. 
In  the  NIR  region,  I  observed  the  zero  phonon  line  with  a  small 
triplet  splitting  at  1282.8nm,  1286.5  and  1287. Inm,  possibly  due 
to  the  spin-orbit  interaction.  The  corresponding  phonon  sideband 
shows  that  at  least  one  phonon  energy  of  107cm*  ^  is  involved.  The 
same  characteristics  also  appear  in  the  peaks  due  to  the  assigned 
transition  3a2->3Ti  of  Cr^+  in  the  Ge  site.  The  average  phonon 
energies  were  found  to  be  644cm-l  and  228cm*  ^  for  the 
transitions  with  energies  15723cm**  and  13736cm**. 

With  the  assumption  of  linear  coupling  between  electronic 
states  and  lattice  vibrations,  we  can  use  the  quantum-mechanical 
Franck-Condon  Principle  to  calculate  the  intensities  of  the  zero 
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Fig.  4.13  Room  temperature  polarized  absorption  spectra 
of  CnBaiMgGeiO?. 
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Fig.  4.14  lOK  polarized  absorption  spectra 
CrrBaiMgGeiOT. 
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phonon  line  and  vibronic  transitions. ^5. 16. 17  According  to  the 
Huang-Rhys  calculation,  the  relative  integrated  intensity  of  zero 
phonon  line  at  T=0  is  expressed  by  the  Huang-Rhys  factor  So 

(Itot/Izpl)T=0=  exp(So)  (4.1) 

where,  Izpi  and  Itot  are  the  corresponding  integral  intensities  of 
zero  phonon  line  and  phonon  sideband,  and  So==£sPs/Ao>s  is  the 
sum  of  Stokes  losses  expressed  in  phonon  numbers;  Ps  is  the 
Stokes  shift  in  energy  for  the  s-th  mode;  2Ps=msO)s^<lso^.  where  mg 
is  mass,  cog  is  the  frequency  and  qos  is  the  shift  of  the  equilibrium 
position  for  the  oscillator. 

As  the  temperature  is  raised,  the  relative  integrated 
intensity  of  zero  phonon  line  is 

(Itot/Izpl)T=  exp(ST)=exp(IsPgCOth(ftci)s/2kT)M(Og)  (4.2) 

By  using  the  above  formula,  (Itot/Izpl)T*lOK  for  the  transition  near 
1287nm  was  found  to  be  58.285.  This  value  corresponds  to 
S(T=10K)=4.  This  big  S  shows  the  large  Stokes  loss  between  and 
3T2  states.  The  Huang-Rhys  factor  of  the  3a2'^^Ti  transition  at 
636nm  was  found  to  be  1.7  at  lOK. 

By  using  the  low  temperature  absorption  data,  the  cubic- 
field  approximation  parameters  are  readily  obtained  by  solving 
the  eigenenergy  equations  in  Appendix  to  be: 

Dq=-91 4.5cm*  1  ,  B=498cm-l  ,  C==4.4B 

The  calculated  and  observed  energy  levels  are  shown  in  Fig.  4.15. 
The  agreement  is  reasonably  good. 

No  emission  was  found  in  Cr:Ca2MgSi207  at  either  room 
temperature  nor  liquid  helium  temperature.  Based  on  the 
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Fig.  4.15  Calculated  and  observed  energy  levels  of  Cr^ 
ions  in  BazMgGeiO?. 

sensitivity  of  the  spectrometer,  the  Cr;Ca2MgSi20  7  emission 
intensity  must  be  at  least  a  factor  of  100  less  than  that  of 
Criforsterite 


4.2.4  Comparison 

Here,  some  results  are  summarized  from  section  4.2.1 -4.2.3. 
First,  Table  4.1  shows  the  absorption  peaks  of  three  Cr  doped 
crystals. 

We  notice  that  all  the  peaks  fall  in  the  same  frequency 
region  especially  for  Cr:forsterite  and  Cr:Ca2MgSi20 ?•  The 
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Table  4.1  Absorption  peaks  of  three  Cr-doped  crystals. 


Crystals 

Absorption  Peaks  (nm) 

Cr;Mg2Si04 

571,  662,  740 

Cr:Ca2MgSi207 

579,  660,  767 

Cr:Ba2MgGe207 

636,  728,  832 

absorption  peaks  of  CT;Ba2MgGe207  are  red-shifted  relative  to  the 
other  two  crystals.  This  is  due  to  the  smaller  energy  splitting 
resulting  from  the  smaller  crystal  field  strength  of  the  larger  Ge 
tetrahedral  site  of  Ba2MgGe207.  These  absorption  characteristics 
must  come  from  having  similar  species  in  nearly  identical  local 
environments  in  these  crystals.  Comparing  the  crystal  structure  of 
forsterite  with  ^kermanite,  the  only  similarity  is  the  Si(Ge) 
tetrahedral  site  with  symmetry  Cs-  Also,  in  the  ikermanite 
crystals  ,  there  are  no  chromium  charge  compensating  species; 
thus,  only  Cr^'*’  can  substitute  for  the  Si^+  or  Ge'^*  sites.  Therefore, 
the  assumption  that  tetrahedrally  coordinated  Cr^'*'  ions  are  in  all 
three  crystals,  Mg2Si04:Cr,  Ba2MgGe207:Cr  and  Ca2MgSi207:Cr  is 
consistent  with  the  room  temperature  absorption  data. 

However,  considerable  differences  in  the  low  temperature 
NIR  absorption  spectra  exist.  First,  zero-phonon  lines  are  absent  in 
Cr:Ca2MgSi207  for  the  transition  3a2->^T2.  Second,  although  zero- 
phonon  lines  were  observed  both  in  the  Cr:Ba2MgGe207(1287nm) 
and  Cr;forsterite(1090nm),  they  occur  in  different  polarizations 
since  the  c-axis  of  dkermanite  corresponds  to  a-axis  of  forsterite 
when  considering  the  orientation  for  the  tetrahedra  (referring  Fig. 
2.4).  The  Crrforsterite’s  relative  absorption  of  the  zero-phonon 
line  to  the  vibronics  is  considerable  larger  than  that  of 
Cr:Ba2MgGe207.  And  most  important  of  all,  no  emission  has  been 
found  from  either  chromium  doped  4kermanite  crystals. 


These  results  did  not  prove  that  the  transition  of  Cr^*  ions  in 
the  Si  site  of  forsterite  is  the  mechanism  for  the  NIR  laser  action 
of  Criforsterite.  However,  the  excitation  spectra  of  Cr:Mg2Si04  at 
20  K  while  monitoring  emission  in  the  NIR  region  traced  almost 
exactly  the  large  absorption  peaks  in  the  visible  region.  But  the 
chances  are  that  even  at  20K,  the  energy  transfer  between  species 
is  not  diminished.  Then,  another  evidence  by  Ref.  14  shows  that 
the  absorption  coefficients  in  the  NIR  and  visible  region  have  the 
same  increasing/decreasing  trends  when  the  forsterite  crystal  was 
annealed  or  oxidized.  This  can  definitely  relate  the  two  spectral 
regions,  or  at  the  very  least,  we  can  say  that  the  laser  transition  of 
forsterite  is  related  to  Cr^'*'. 

Let  us  compare  crystal  structures  between  forsterite  and 
Akermanite.  The  local  environments  are  so  similar  in  that  they 
both  have  the  Cs  symmetry  and  a  Si-0  average  bond  length  1.63 A. 
But  we  must  also  consider  the  whole  structures  of  the  crystals 
because  the  d  electron  wavefunctions  have  strong  interaction  with 
the  lattice  vibrations.  In  considering  the  whole  crystal  structures, 
forsterite  and  Akermanite  have  some  differences. 

We  can  use  the  simple  configuration  coordinate  model  to 
compare  the  electron -lattice  coupling  of  the  three  crystals. 

By  using  Eqs.  4.1  and  4.2,  the  corresponding  Huang-Rhys 
factors  for  Cr .-forsterite  and  Cr:Ca2MgSi207  can  be  obtained.  For  Cr: 
forsterite,  the  NIR  transition  at  I.lp.m  has  Izpi/ltoi~  0.2927  and 
S=1.2286  at  lOK  as  compared  to  Cr:Ba2MgGe207  (S=4)  and 
Cr:Ca2MgSi20 7(-several  tens)  of  some  transitions.  So  for 
Cr:Ba2MgGe207  Cr:Ca2MgSi207,  the  coupling  is  strong(S»l)^5  and 
the  equilibrium  position  shift(Stokes  shift)  between  ground  state 
and  excited  state  is  large.  As  shown  in  Fig.  4.16(a),  when  the 
Stokes  shift  is  large,  the  overlapping  of  the  phonon  wavefunctions 
between  the  ground  and  excited  state  is  small,  as  is  the  zero- 
phonon  line  intensity. 
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If  the  Stokes  shift  is  so  large  that  one  branch  of  the  ground 
state  potential  energy  parabola  intersects  the  excited  state 
potential  parabola,  this  intersection  will  allow  a  non-radiative 
relaxation  path  shown  in  Fig.  4.16(b).  Therefore,  non-radiative 
relaxation  will  dominate  and  no  radiative  relaxation  can  be 
detected.  This  may  be  the  reason  why  Cr:Ca2MgSi20 7  and 
Cr:Ba2MgGe207  have  none  or  very  weak  emission. 
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4.3  Results  and  Analysis  of  EPR  Spectroscopy 
4.3.1  Cr:Akermanite(Ba2MgGe207)^ 

In  EPR  spectroscopy,  0.5%  and  0.01%  chromium  doped 
Ba2MgGe20  7  were  examined.  Three  different  signals  were 
detected.  Two  of  those  signals  which  do  not  scale  according  to  the 
chromium  concentration  are  from  Mn2+  and  Fe3+  impurities. 

The  one  that  scaled  according  to  the  Cr  concentration  must 
be  due  to  chromium  ions.  With  the  magnetic  field  in  an  arbitrary 
direction,  we  can  observe  four  EPR  signals  corresponding  to  four 
crystallographically  equivalent  magnetic  sites.  These  four  sites  are 
shown  in  Fig.  4.17. 

When  magnetic  field  points  to  a  specific  direction,  some  of 
these  sites  will  become  magnetically  equivalent  in  a  manner  that 
indicates  the  incorporation  of  this  species  at  the  crystallographic 
4e  sites  of  Ba2MgGe207.  Both  barium  and  germanium  lie  at  4e 
sites.  But  due  to  the  large  disparity  of  barium  and  chromium  ionic 
radii(0.44A  vs  1.42A),  it  is  unlikely  that  chromium  ions 
incorporate  into  Ba  site.  So  it  is  inferred  that  chromium  ions 
occupy  the  Ge's  tetrahedral  sites. 
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Fig.  4.16  (a)  Configuration*coordinate  diagram  showing 
the  zero'phonon  line  transition,  (b)  Configuration* 
coordinate  diagram  for  non*radiative  processes. 


Fig.  4.17  BaiMgGeiO?  P4  2im  space  group  symmetry. 
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The  field  for  resonance  curve  which  is  plotted  along  several 
directions  is  shown  in  Fig.  4.18.  This  spectral  feature  is  attributed 
to  Cr^'*’  in  Ge  site. 

Only  one  of  three  possible  microwave  transitions  was 
detected  indicating  large  zero-field  splitting.  The  observed 
transition  is  the  one  from  the  Sz~-1  state  to  the  Sz=:l  state  from 
Fig.  3.  3.  The  angular  variation  of  the  field  for  resonance  curve 
was  well  described  by  the  usual  spin-hamiltonian  for  a  spin  one 
system  described  in  Section  3.3.1. 

rt  =  pH.g.S+D[Sz*-l/3S2]+E(Sx*-Sy*)  (4.3) 

In  Section  3.3.1,  the  above  formula  is  used  for  a  tetragonal  crystal, 
but  it  can  also  be  applied  to  any  crystal  with  lower  symmetry  (in 
this  case,  Cs  symmetry  for  Cr^"*"  in  the  Si(Ge)  sites  of  forsterite  and 
4kermanite),  because  we  can  always  diagonalize  the  spin- 
hamiltonian  by  eliminating  the  SxSy+SySx*  SxSz+SzSx  and 
SySz+SzSy  terms.  After  diagonalization,  the  axes  of  the  spin- 
hamiltonian  may  not  be  parallel  with  the  crystal  axes.  For  Cs 
symmetry,  the  only  symmetry  operation  is  a  mirror  plane.  One  of 
the  spin-matrix  axes  must  be  perpendicular  to  the  mirror  plane, 
and  the  other  two  lie  on  the  mirror  plane  with  unknown 
orientations. 

After  fitting  the  spectra  to  spin  hamiltonian  Eqs.  4.1  with 
the  least  square  method,  the  fine  structure  parameters  values  are 
(standard  deviations  in  parentheses): 

101=1.1X102(0.3X102)  GHz,  IEI=4.1(0.1)GHz 

The  dominant  fine-structure  axis  (z  axis)  lies  in  the  mirror 
plane(a±b  plane)  and  makes  an  angle  of  approximate  7r/3  with  the 
crystallographic  c  axis,  while  the  minor  (x  axis)  is  normal  to  the 
mirror  plane. 
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Fig.  4.18  Angular  dependence  of  resonant  EPR  field  in 
doped  BaiMgGeiO?.  (After  Ref.  18) 
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Magnetic  Field  Orientation 


4.3.2  Cr:Forsteritei8 


The  EPR  spectra  of  Cr:forsteriie  are  more  complex  than 
Cr:Ba2MgGe207  due  to  the  fact  that  charge  compensated  chromium 
ions  of  various  valences  as  well  as  other  impurities  are  present  in 
Cr'.forsterite.  The  dominant  signal  which  is  due  to  Cr3+  in  Mg  site  is 
readily  identified^  ^ .  With  the  assistance  of  the  known 
Cr:Ba2MgGe207  system,  one  signal  very  similar  to  the  Cr^*  center 
in  Ba2MgGe207  was  detected.  The  field  for  resonance  as  a  function 
of  magnetic  orientation  is  shown  in  Fig.  4.19.  This  signal  pattern  is 
consistent  with  the  4c  site  in  forsterite  and  is  thus  consistent  with 
Cr4+  substitution  for  silicon  (referring  the  Table  2.1). 

After  fitting  to  Eqs.  4.3,  the  fine  structure  parameters  are 
obtained: 

lDI=61(l)GHz  and  1E1=4.588(0.003)GH2 

The  dominant  axis  (z)  lies  in  the  mirror  plane  and  makes  an  angle 
of  rt/4  with  crystallographic  a  axis  and  the  y  axis  is  normal  to  the 
mirror  plane. 


4.3.3  Comparison 

From  EPR  spectra,  Cr4+  ions  in  tetrahedral  sites  are  identified 
in  both  Cr:Ba2MgGe207  and  Cr'.forsterite.  These  results  are 
consistent  with  the  results  obtained  from  the  optical  methods. 
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Fig.  4.19  Angular  dependence  of  resonant  EPR 
doped  forsterite.  (After  Ref.  18) 
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Chapter  5 

Final  Discussions 


5.1  Summary 

(1)  Czochralski  growth  of  chromium  doped  Ca2MgSi207  and 
Ba2MgGe207  were  developed.  Crystals  were  grown  either  by  rf 
induction  heated  or  resistance  heated  furnaces.  The  different 
aspects  of  three  furnaces  were  explored. 

(2)  Chromium  doped  ikermanites  (Ca2MgSi207  and  Ba2MgGe207) 
were  grown  to  assist  the  identification  of  the  impurity  centers  in 
chromium  doped  forsterite  crystals.  Large  absorption  peaks  in  the 
visible  frequency  regions  are  easily  identified  in  all  the  samples. 
This  similarity  and  the  unique  environment  of  the  dkermanite 
crystals  which  allow  only  Cr^'*'  going  into  the  Si(Ge)  site  confirm 
the  assumption  of  Cr^'*'  ions  in  forsterite. 

(3)  Electron  paramagnetic  resonance  spectroscopy  directly 
measure  the  ground  state  splitting  and  site  symmetry  of  the 
impurity  centers.  The  results  are  consistent  with  the  assumption 
of  Cr*"**  substituting  for  the  Si(Ge)  tetrahedral  sites. 

(4)  From  both  optical  and  magnetic  resonance  methods,  only  Cr^'*' 
ions  in  the  silicon(germanium)  tetrahedral  sites  have  been  found 
in  chromium  doped  &kermanite.  Chromium  ions  in  forsterite  not 
only  occupy  the  Mg(Ci)  and  Mg(Cs)  octahedral  sites,  but  also 
substitute  for  the  silicon  tetrahedral  sites. 
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(5)  There  exist  considerable  differences  between  the  low 
temperature  absorption  spectra  of  Crrforsterite  and  Cr:4kermanite. 
Additionally  no  emission  was  detected  in  the  chromium  doped 
ikermanite  samples.  This  leads  us  to  the  conclusion  that  the 
placement  of  Cr^+  in  tetrahedral  sites  does  not  guarantee  a  laser 
transition  in  the  near-infrared  region. 

(6)  A  possible  explanation  why  Cr:4kermanite  do  not  have 
emission  is  given,  although  other  possibilities  are  not  ruled  out. 
For  example,  the  coupling  between  Cr^-*-  and  Cr*+  or  interaction 
between  charge  compensation  ions  are  also  possible.  But  at  least, 
the  laser  action  in  Criforsterite  is  Cr4+  related. 


S.2  Survey  of  Configuration  in  Tetrahedral  Sites 

During  the  period  of  this  work,  several  other  potential  solid 
state  laser  materials  using  the  d^  configuration  in  tetrahedral  sites 
were  developed.  Those  materials,  along  with  Cr:forsterite  and 
Cr:4kermanite  are  discussed  in  the  following  paragraphs: 

(1)  Cr:Mg2Si04:  Cr:forsterite  have  both  Cr^-*-  in  octahedral  sites  and 

in  tetrahedral  sites.  Cr2+  ions  may  also  substitute  for  Mg2+. 
Laser  tuning  ranges  from  1170nm  to  1330nm.  Both  continuous 
and  pulse  operations  have  been  demonstrated.  Mode-locked  and 
Q-switched  lasing  mode  have  also  been  achieved. 

(2)  Cr  doped  melilite  group:  several  crystals  tested  include  the 
Akermanites  (Ca2MgSi207.  Ba2MgGe207)  used  in  this  study, 
g eh  leni  te ( C a2  A  1 2  S  i  O  7,  C  a  2G  a  2 S  i  O  7 )^  .  and  mixed 
systems(Ca2(MgxAlySii.x-y)307).  Gehlenites  have  the  same 
structure  as  ikermanites  with  aluminum  replacing  magnesium 
and  half  of  the  silicon  atoms.  It  was  found  that  the  intensity  of 
luminescence  of  chromium  ions  in  the  near  IR-region  significantly 
decreased  with  the  change  of  host  composition  from  aluminium 
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content  melilite  to  magnesium-content  melilite.3  Although  the 
emission  spectra  observed  in  Cr:gehlenite  are  likely  due  to  Cr^  + 
ions,  an  emission  band  extending  below  lOOOnm  (in  Cr:Ca2Ga2Si07) 
and  a  nonexponential  fluorescence  decay  (in  both  Cr:Ca2Al2Si07 
and  Cr;Ca2Ga2Si07)  can  not  be  explained  by  solely  Cr^'*'  ions  in 
tetrahedral  sites.  No  laser  action  has  been  reported  yet. 

(3)  Chromium  doped  garnets3.4,5,6.7,8.9;  garnet  crystals  have  space 
group  symmetry  Ia3d.  When  chromium  ions  were  doped  into 
garnets,  they  are  likely  to  enter  into  both  octahedral  and 
tetrahedral  sites.  Tetrahedral  sites  in  garnets  have  the  S4  cubic 
symmetry  which  is  reflected  in  the  absorption  spectra.  The  same 
characteristic  3a2->3Ti  absorption  peaks  are  also  found  in  the 
absorption  spectra.  Chromium  doped  yttrium  aluminum  garnet 
(YAG)  have  been  demonstrated  to  lase  between  1.3pm-1.5|im.  The 
lack  of  zero-phonon  line  is  due  to  the  forbidden  transition 
between  3A2  and  3T2  levels.  The  fluorescence  :*ecay  is  exponential 
with  a  lifetime  of  3.6  ps.  Similar  emission  s^  •  were  also  found 
in  Cr:Gd3Sc2Ga30 12  and  Cr:Y3Sc2Ga30i2.  l-t  none  of  these  were 
demonstrated  to  lase,  probably  due  to  the  excited  state  absorption 
3T2->3Ti. 


(4)  Mn5+  doped  Sr5(P04)3Cl  and  Ca2P04Cn®:  Sr5(P04)3Cl  has  space 
group  symmetry  P63/m,  and  Ca2P04Cl  symmetry  Pbcm.  The 
tetrahedral  sites  have  respectively  C2  and  Cs  symmetry  in 
Sr5(P04)3Cl  and  Ca2P04Cl.  The  Mn5+  ion  is  most  likely  stabilized  in 
the  tetraoxo  coordination  by  isomorphous  substitution  for  the  P5  + 
ion,  giving  MnOa^'  tetrahedra.  The  room  temperature  absorption 
spectra  show  large  peaks  which  indicates  that  strong  coupling 
between  local  vibrations  and  electronic  transitions  probably 
exists.  The  coupling  to  the  local  vibrations  would  be  no  surprise 
because  the  high  valence  state  (5+)  of  manganese.  Fluorescence 
decay  is  double  exponential  and  has  much  longer  lifetimes 
compared  to  that  of  Cr:Mg2Si04  and  CrrYAG. 
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The  similar  behavior  of  probable  Mn5  +  (d2)  in  tetrahedral 
sites  should  give  us  some  confidence  about  the  assumption  of 
laser  transition  of  Cr^+(d2)  jn  tetrahedral  sites.  However  a  lot  of 
questions  still  need  to  be  answered.  It  seems  that  emission  does 
occur  when  there  are  many  different  substitutional  sites  available 
for  chromium.  So  the  possibility  of  cross  interaction  can  not  be 
ignored. 


5.3  Future  Research 

(1)  A  possible  experiment  by  combining  optical  excitation  and 
electronic  paramagnetic  resonance  can  be  used  to  directly 
measure  properties  of  the  excited  state  of  the  laser  transition  at 
l.l-1.3|im. 


(2)  Observe  the  shift  of  the  zero  phonon  line  in  the  near  infrared 
region  when  applying  external  magnetic  and  electric  field.  This 
should  give  us  more  information  about  the  site  symmetry  and  the 
energy  levels. 

(3)  Grow  other  kinds  of  crystals  with  chromium  in  tetrahedral 
sites  to  study  the  coupling  effect  between  lattice  vibrations  and 
electronic  transitions  in  different  environments. 
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Appendix 

Matrix  Elements  and 
Eigenenergies  of  d2  and  d3 
configurations 


The  parameters  used  in  the  following  sections  are  defined  as 
followed. 

Racah  parameters: 

B=:F2-5F4 

C=35F4 

Slater-Condon  parameters: 

F2=F2/49,  F4=F4/441 

Flf(dd)=/ri2dri/r2^dr2Rd2(ri)R(i2(r2)r<*^/r>*^+^ 

where  r<  is  the  lesser  and  r>  is  the  greater  of  n  and  r2,  Rd(ri)  and 
Rd(r2)  are  the  radial  functions  of  3d  electrons. 

When  calculating  the  eigenenergies,  the  appropriate  Dq  values 
must  be  added  to  diagonal  terms  of  the  matrices. 


A.l  d2  configurations 


A.1.1  Matrix  Elements 
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A. 1.2  Eigenenergies(E/B) 

lAi  (1/2)  {18+9  C/B±[I00+25  (C/B)2+400  (Dq/B)2+100  C/B- 
80Dq/B-40  (C/B)  (Dq/B)Jl/2} 

3A2  -8+10  Dq/B 

lE  (1/2)  {1+4  C/B±[49+400  (Dq/B)2.40  Dq/B]l/2} 

iTi  4+2  C/B 

3Ti  (1/2)  {-1-10  Dq/B±[225+100  (Dq/B)2+180  Dq/BlI/2} 

'T2  (1/2)  {1+4  C/B-10  Dq/B±[49+100  (Dq/B )2.20Dq/Bl  1/2) 

3T2  -  8 


A.2  Configuration 
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A.2.1  Matrix  Elements 
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A. 2. 2  Some  of  Eigenenergies(E/B) 
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